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ASSUMING LIGHT TRACES MASS: LENSPERFECT ANALYSIS OF ABELL 1689 
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Draft 

ABSTRACT 

We present a strong lensing mass model of Abell 1689 which resolves substructures ~ 25 kpc across 
(including about ten individual galaxy subhalos) within the central ~ 400 kpc diameter. We achieve 
this resolution by perfectly reproducing the observed (strongly lensed) input positions of 168 multiple 
images of 55 knots residing within 135 images of 42 galaxies. Our model makes no assumptions about 
light tracing mass, yet we reproduce the brightest visible structures with some slight deviations. A1689 
remains one of the strongest known lenses on the sky, with an Einstein radius of Re — 47.0" ± 1.2" 
(143^4 kpc) for a lensed source at Zg = 2. We find a single NEW or Sersic profile yields a good fit 



simultaneously (with only slight tension) to both our strong lensing (SL) mass model and published 
weak lensing (WL) measurements at larger radius (out to the virial radius). According to this NEW 

^°J X 10i^Mf^/i7n^ (Mann = 1.8l:n1 x 10^^ Mo,h:;n) within the virial 



fit, A1689 has a mass of M„,^ = 2.0^^1 x lO^^Mm/i 



o"7n 



11.5 



+ 1.5 
-1.4 



radius r^ir = 3.0±0.2 Mpc hj^ {r2oo = 2.4itn'2 Mpc hj^), and a central concentration c„ 
(c2no = 9.2 ± 1.2). Our SL model prefers slightly higher concentrations than previous SL models, 
bringing our SL+WL constraints in line with other recent derivations. Our results support those of 
previous studies which find A1689 has either an anomalously large concentration or significant extra 
mass along the line of sight (perhaps in part due to triaxiality) . If clusters are generally found to have 
higher concentrations than realized in simulations, this could indicate they formed earlier, perhaps as 
a result of early dark energy. 

Subject headings: gravitational lensing — methods: data analysis — galaxies: clusters: general — 
galaxies: clusters: individual (Abell 1689) — cosmology: dark matter 



1. INTRODUCTION 

Galaxy clusters are the largest gravitationally bound 
objects yet to form in our universe and prove interesting 
to study both intrinsically and toward other ends. Maps 
of their dark ma tter and bary ons yiel d insights into struc- 



controversial) findings, from the Bullet Cluster (Marke 



vitch et al.|2004[ [Clowe et al. 2006| 
the "dark matter ring" in CLU024 



Randall et al. 

ee ct al. 2007^ 

(Mahdavi et al 



to 
and 



the " cosmic train wreck" in Abell 

2007|). These results are exciting in part because they 



ture formation (Umetsu et al. 2009; Kawaharada et al 



|2010| and can even constrain the nature of dark matter 
particles (specifical l y, their self-inter acting cross-section; 



Randall et al.l2008[ [Fmg et al.|2009| . Cosmological con 
straints can be obtained troni nurnber counts of clus 



ters with measured masses dMantz et al. 2009; Vikhlinin 



et al.|[2009 Rozo et aL]|2010p or from the observed grav- 



itational l ensing of more dis tant objects with measured 
redshifts (Gilmore & Natarajan 2009). As the strongest 
gravitational lenses on our sky, clusters can also be used 
as cosmic telescopes allowing us to peer further back in 
time to reveal galaxies in the dista nt ( z > 7) universe 
( [Bradley et"al| 2008 ; Ri chard ct al. 2d08| or study galax- 
ies at "modest redshifts [z ^ 3 to 4) in greater detail 



( [Bunker et aT][2000, ,Frye et al.[[2007 l. 

By mapping the dark matter distributions within 
galaxy clusters (as projected on our sky), gravitational 
lensing analyses have yielded many exciting (and at times 
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reveal the distribution of dark matter without assuming 
that light traces mass (hereafter, "LTM'rJ). 

Similarly exciting dark matter maps can now be ob- 
tained in finer detail in galaxy cluster cores thanks to 
strong gravitational lensing (hereafter, SL) analysis of 
high quality data. The galaxy cluster Abell 1689 is one 
of the strongest gravitational lenses on our skjF] with 
an Einstein radius of Re ^ 47" for a backgrouna object 
at Zs = 2. Deep (20-orbit) multiband ACS GTO ob- 
servations of this cluster reveal over 100 str ongly-lensed 
multiple ima ges of 30 background galaxies (Broadhurst 
et al. 2005a), a huge leap forward in the number of 
bL constraints available for any single massive body. 
Based on these multiple images (and others identified 
since), we derive a mass model which resolves substruc- 
ture approximately 25 kpc (8") across within the central 
^ 400 kpc (2') diameter. This is the highest resolution 
mass model to date of any galaxy cluster without as- 
suming LTM. For comparison, weak lensing (hereafter, 
WL) mass maps typically resolve structure on ~ 1' scales 
from ground-based imaging or ~ 45" from space ( Hey-| 



We will avoid describing models as "paramctric]]or|^non- 



parametrTc"7^ince''s^rTctly speaking all mass models SoTTave'^er 



rameters. H ere we introduce the more importa nt distinction be 
t ween "L'lTl" and "non-LTM" models. ^ 

'^ !riee|/^itrm et al.|l|2009a[| for the current strongest gravitational 



lens. 



Coe et al. 



.|2008[ ) 
3Des cIi 



mans et al.l2008|. Of course WL complements SL nicely, 



as WL probes cluster mass distributions to much larger 
radii (the virial radius r^r > 2 Mpc) where the lensing 
strength is weaker. 

Traditional SL analysis methods were not able to fully 
process the large numbers of multiple images revealed in 
the ACS images of A168 9. LTM models produced by 



Broad hurst et al 



highly magnified background galaxies and multiple im- 
age systems. 

The observations (obtained in 2002 June) delivered, 
and the 100+ multiple images still far surpass any other 
gravitational lens. To date, other clusters observed to the 
same depth have yielded far fewer multiple images: 53, 
35, and 33 froin A1703, A2218, and CL0024 respectively 



J2005al h ereafter B05), [Z ekser et al., ,, . 

he reafter ^06), I Halkola et al. ( 2U06 hereafter 2009b|). A1689 s many SL constraints allow for detailed 



H06), and |Limousin et al.| ( |2007[ hereafter LOt) all failed 



( Richard etaLj2009a[|Eliasd6ttir et al.|2007[|Zitrin et al. 
z009bk. A1689 s many SL constraints allow for detailed 



to reproduce the 100+ observed multiple image positions 
by > 2.5" (see Table [l|. This is roughly 50 times the ob- 
servational uncertainties of one pixel or so (~ 0.05"). 
By failing to match the tight observational constraints, 
these models discard a great deal of information available 
in the quality ACS images. 

Despite their greater flexibility, non-LTM models were 
unable to reduce these resid uals given all 100+ multiple 



images. 

uals similar to those in Ll'M models. iSaHa et al. 



Diego et al. ( 2005b I, using SL AH^ allowed resid 



5 ana 



(2006), 



using PixeLens, produced non-LTM models which per 
fectly reproduce some of the data, but computational 
limitations restrict PixeLen s to fitting only 30 m ultiple 
images at a time. Similarly, Julio & Kneib (2009), using 
a hybrid scheme combining L'i'M and non-Ll M substruc- 
ture, leave relatively small (0.28") average residuals, but 
only for a subset of 28 multiple images. 

In this paper, we present non-LTM mass models which 
perfectljrl reproduce the observed (lensed) input posi- 
tions of 1d8 multiple images of 55 knots residing within 
135 images of 42 galaxies. The resolution of a recon- 
structed mass model is given roughly b y the spatial den- 



sity of th e multiple image constraints ( Coe et al. |2008 
Coe|2009 ) . Thus our model which incorjaorates 135 mul 
tiple images has ~ 4x greater spatial resolution (^ 2x 
greater along each axis) than the aforementioned models 
which incorporate only ~ 30 multiple i mages. 

Our SL analysis method (LensPerfect, Coe et al.|2008 l 
was made possible just recently thanks to the develop- 
ment of mathematical tools enabling one to obtain curl- 
free interpolations of a vector field which perf ectly re- 
produce th e input data at the input positions (Fuselier 



2006[ |2007[ ) . We make only minimal assumptions about 
the physicality of the mass distribution, including no as- 
sumptions about LTM. 

As the strongest known gravitational lens on our sky 
for some time, A1689 has provided us with an excellent 
laboratory. It was the first cluster to have i ts mass mea- 
sured via both analyses of weak lensing (Tyson et al. 
1990 1 a nd magnification of the background galaxy pop 



ulation (Taylor et al. 1998). Yet before A1689 was ob- 



served with ACS, no multiple images had been identified 
in the field except for two fold a rcs, our 8ab and 13abc 
( Miralda-Escude & Babul 1995). These arcs were suf- 
hcient for a rudimentary strong lensing analysis to be 
performed yielding an Einstein radius of Re ~ 45", the 
largest of any known lens at the time. With such a large 
Einstein radius, the ACS GTO team devoted 20 orbits to 
imaging the cluster, confident that it would reveal many 

"^ Strong Lensing Analysis Package 

® While the observed image positions do have (small) obser- 
vational uncertainties of a pixel (0.05") or so, our solutions do 
perfectly fit the data as input. 



mass modeling, spawning many publications and helping 
to make A1689 one of the best studied galaxy clusters. 

One result that stands out from both lensing and X-ray 
analyses of A1689 is that its mass appears to be more cen- 
trally concentrated than predicted by CDM simulations 
of structure formation. Simulated dark matter halos have 
mass profiles which a re generally well desc ribed by NEW 



([Navarro et al.||1996^ or ^_ 

2004[) profiles (e.g., Navarro et al.||2010|). More massive 



Einasto / Sersic (Navarro et al. 



clusters (such as A1B89) finished forming later when the 
universe was less dense overall, and thus are found (at 
least in simulations) to be le ss centrally conce ntrated. 
Yet where A1689 is expected ( |Duffy et all2008| to have 
a NEW concentration of C200 = 3^0T^TTcvir = ^-Qti^), 
it has been shown observationally (Table [2]) to have a 
much higher concentration C200 ^7 — 11 {cyir '-^ 8— 12)rl 
As one of the strongest lenses on our sky, we might ex- 
pect Al689's concentration to be on t he high side. How- 



2005 



ever, even acco unting for triaxiality (lOg uri et al 
Corless et al.|200 9 ) and selection (lensing) bias dhLennawi 
[ Oguri fc Bla ndford'^ 2009[ i Meneghetti et al. 



et al. 



2010|, A 1689's high conc entrati on stil l seems extremely 



unlikely (Broa dhurst fc Barkana| 2008 1 . 

Lensing analyses of larger samples of clusters seem 
to support t he idea that clusters may have higher con- 
centrations CBroad hurst et al ||2008[ (Oguri et al.'l2009 
ISereno eTal . 2010) and larger Einstein rad ii (Broadhur^ 
fe Barkana(2008[ [Richard et al.|200"9b{|Zitrin et al |2010a| 



3tl 



than simulated clusters. However only a small lensing- 
biased sample of these have been studied in sufficient 
detail. We note that joint SL+WL fitting is required to 
constrain mass profiles we ll and concentrations to ; ^ 10% 
according to simulations (Meneghetti et al. 2009). The 
addition of X-ray, SZ, and/or velocity disp ersion data 



can constrain the m ass profiles f urther still (i Sand et al. 
Newman et al.,,2009, Lcmze"et al.| |2009fT 



2008 



More conclusive results are expected from the 
"CLASH'p] Treasury Project, a large (524 orbit) HST 
Multi- Cycle Treasury program (P.I. Postman) to observe 
25 X-ray-selected galaxy clusters at 0.18 < z < 0.9, each 
to a depth of 20 orbits, or the depth of the ACS GTO 
images of five clusters including A1689. Combined with 
Subaru images and other data, these observations should 
yield conclusive results (see Q. 

Might baryons, lacking from the simulations discussed 
above, be responsible for higher mass concentrations in 



^ [Morandi et al.| | |2010| claim to account for A1689's high con- 
centration (along with the discrepancy between lensing and X-ray 
masses) by fitting an 3-D ellipsoidal gNFW profile (with vari- 
able central slope) simultaneously to the SL, WL, and X-ray data. 
While the method is impressive, simulated clusters have yet to be 
analyzed in the same way (they are normally fit to spherical NFW 
profiles), so it is unclear that a direct comparison can be made. 

^^ C luster Lensing And Supernova survey with Hubble, 
|http://www. stsci.edu/~postman/CLASH/ 
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TABLE 1 
Published strong lensing mass reconstructions of A1689 based on the ACS images 



Alias 

B05 

Z06 

H06 
HOT 
L07 

LP 



Paper 

Broadhurst ct al. (2005al 
DiCKO ct al. ( 2UU5b ) 
Zekscr ct al. ( ^UUCT 
Saha ct ai. ( 2uub1r 
Halkola ct al. {2{M6] 
Halkola ct al. ( 2UUY ) 
Limousin ct al. (2UU7|| 
Jiillo & Kncib (2009) 
T^is paper 









Average 


residuals 


Assume LTM 


Galaxies 


Multiple Images 


source plane 


image plane 


/ 


30 


106 




3.2" 


X 


30 


106 




3" 


/ 


22-30 


74-106 


0.7"-1.2" 
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7,9 


28,30 








/ 


31+arc 


107 




2.7" 


/ 


31+arc 


107 




2.5" 


/ 


34 


113 


0.45" 


2.87" 


/x^ 


12 


28 




0.28" 


X 


42 


135 









'^ Observational uncert ainties are ~ 0.05" in the image plane. 

Julio &; Kneib (20091 include both LTM and non-LTM substructure components. 



TABLE 2 
Published concentration measurements of A1689 



4.7 
5.7 
7.9 

q c+O.S 

5.2 ±0.3 

°-^-1.6 

D+1.2 

'-0.8 

7+0.34 

-0.5 

6 ±0.5 



lO.i 

5.7! 



10. 7"* 



-g+0.3 
'•"-0.5 

6.0±0.6<= 
7.6 ±1.6 
7.6 ± 1.3 
9.2 ± 1.2 

7+4.5 
-2.7 

10 1+"-* 

^'^■^-0.7 
12.5t|? 

9 8+°-^' 
'^•'^-0.6 

11.10 
12.2 ± 6. 7^* 

7+1.7 



2.6 



+ 1.3 



7.7] 

5.6 

5.3_j 2 
6.6 ±0.4 
10.9 ±3.5 

Q 7+0-7 

> 10.4 
4.58 ± 0.34'' 



Alias 



B05 

Z06 
H06 
H06 
L07 
L07 
LP 



Paper 



Clowe fc Schneider | |200T) 



King ct al. (20U2b,) 
King ct al. (20U2ai 
Clowe (2U03) 



Bardcau ct al.J|2005i 
Bardeau ct ai. fSDUfi 
Broadhurst ct al. (2005a) 
Broadhurst et al., t2005b|) 



Zekser et al.' T2006, 
Halkola ct al. (2006i 
Halkola ct al. ( 2006' i 
Limousin ct al. (2007 1 
Limousin ct ai. i20(TT\ 



inis work 
This work 



Umetsu fc Bro adhurst (2008) 
Umetsu & Broadhurs-t^ (^008J 



Umetsu et al. 



(20091 



Ume tsu ct al.'(^009i 
Corl css et al. ( 2UU9;' 
orlcss et al. ( 2009' 
Andcrsson & Madejski 



rtiem cr-aorcnscn cf^ 
Peiig ct al. (2009) 
Peng ct ai. ( 2Uut)|) 



I 2004 1 



isnra 



tusco-i'cmiano et al 



Lemze et al. 



L.emze et al. 



2008 1 



2Dro 



ivioranai et af 



11217 



]l |2009l ) 



0} 



Constraints'' 

WL 

WL 

WL (infrared) 

WL 

WL 

WL 

SL 

SL+WL 

SL 

SL 

SL+WL 

SL 

WL 

SL 

SL+WL 

WL+counts 

SL+WL+counts 

WL 

SL+WL+counts 

WL 

WL + RE + priors 

X-ray 

X-ray 

X-ray (non-parametric) 

X-ray (parametric) 

X-ray 

SL+WL + X-ray 

SL+WL+counts + X-ray - 

SL+WL + X-ray 



dynamical 



Note. — For previous compilations of concentrations derived for A1689, sec 'Comcrford &| 
[Natarajan (2007 Table Al); Limousin ct al. (2007) Table 4); |Umetsu fc Bro adhurst (2008, Table 
TTT lCoriess et ai.| ( |2009| Table 4). 

" Fits to spherical NFW profiles unless indicated otherwise. Concentrations assuming an 

overdensity of 200 (converted from quoted values if necessary) . 

^ SL = strong lensing; WL = weak lensing; RE = Einstein radius; counts = magnification 

depletion; X-ray; dynamical = cluster galaxy velocity measurements 

'^ 3-cr uncertainties 

'^ Includes uncertainties due to halo triaxiality 

'^ Fit to ellipsoidal gNFW profile (with variable central slope) 



Coe et al. 



nature (e.g., Barkana fc Loeb 2009)? R ecent hydro 



dyna mical simulations (Uutty et al 



2010 Mead et al. 



2010 1 show this is unlikely. Baryons appear to mcrease 
cluster concentrations only modestly at best (~ 10%), 
though they more likely decrease concentrations slightly 
(as found when strong AGN feedback is included in the 
simulations) . 

If real clusters in fact have higher concentrations than 
simulated clusters, this could imply that clusters formed 
earlier in nature than in simulations. One mechanism to 
explain such early growth is a small but non-negligible 
amount of dark e nergy in the early universe, say flnp, '^ 
0.10 at z = 6 dFedeh fc Barteimann] |2007[ [Sadeh fc 
Rephaehl [20081 JFrancis et al.||:^0001 |(irossi fc ^pmgel 
2009p . The adaitional dark energy actually suppresses 
formation of structure, but this means that structures 
must have formed earlier to reach the abundances ob- 
served today. Perhaps high cluster concentr ations along 
with detections of massive halos at z > 1 dEisenhardt 



erar|[2008l I Jee et al.|[2009l [Huang et al.||2U09| | Fapovich 



et al. 



tiona 



2010 jSchwope et al.| 2010 ) are providing observa- 
hints of such early dark energy (EDE). 
We note that semi-analytic modeling of cluster forma- 
tion (in a "standard" ACDM universe) suggests that high 
concentrations (c ~ 10) may be fairly common in nature 
and that cluster mass profiles may di ffer slightly from 
the N EW profiles found in simulations ( Lapi & Cavaliere 
[20091). 



In this paper we concentrate on presenting our revised 
multiple image identifications, our mass model, and mass 
profile fits to SL and WL data, including measurements 
of the mass concentration. In future work we will take 
greater advantage of our method's main strength: the 
ability to map substructure without assuming LTM. 

Our outline is as follows. The observations, object de- 
tections, 12-band photometry, and photometric redshifts 
are described in fpl In E|3l we discuss the observed mul- 
tiple image systems, including our additions and revi- 
sions to those previously identified. In fQwe review our 
LensPerfect method and discuss some minor improve- 
ments we have made to it. Our mass maps are presented 
in ^ and the mass profile is analyzed in ^ We discuss 
substructure and future work in 371 and summarize in ^^ 

We use a concordance cosmology of (O^, ^a, h) — 
(0.3, 0.7, 0.7). In this cosmology , I" ^ 3.11 kpc/ t:7n^ 



at A1689's redshift (measured by Frye et al. 2007) of 
z = 0.187. Though not always explicitly noted, distances 
and masses are given in units of h^^ , where Hq = TOh^o 
km/s/Mpc and h = 0.7/i7o. These values would be 30% 
lower if quoted in units of h~^ instead. 

2. PHOTOMETRY AND PHOTOMETRIC REDSHIFTS 

Multiband observations provide color information 
which is absolutely essential to the identification of mul- 
tiple images. They also allow us to obtain photometric 
redshift estimates where spectroscopic redshifts are un- 
available. Redshifts are essential to the mass model, as 
lensing defiections scale with redshift. 

Abell 1689 is among the best studied galaxy clusters, 
having been imaged in 12 broadband filters from the 
near-ultraviolet to near-infrared. This allows us to obtain 
robust photometric redshifts, as we describe below. By 
carefully modeling and subtracting the light from most of 



the cluster galaxies, we recover many faint objects lost in 
their glare, including demagnified central images which 
allow us to constrain the mass model at small radius. 

Much of the analysis described in this section was per- 
formed for and utilized by B05. 

2.1. Observations 

Abell 1689 (RA=13'*11"30?13, Dec=-01°20'16'.'2 
[J2000]) has been observed in four filters {g'r'i'z') with 
ACS and eight more [UBVRIJHKs) from the ground. 
Details are provided in Table |3l 

The ACS observations are among the deepest to date 
for any galaxy cluster. In 2002 June, 20 orbits of HST 
ACS GTO time were used to obtain deep exposures in the 
ff475: ^625, ^775, and zgso passbauds (4, 4, 5, and 7 orbits, 
respectively). More details about the ACS observations 
can be found in BOS. 

2.2. Galaxy Detection 

Objects were detected in an ACS g' + r' ~\- i' + z' de- 
tection image (with each image normalize d to its back- 
ground RMS) using SExtractor (Bertin & Arnouts 1996). 
Our first detection was conservative (DEBLEND_MINAREA 
= 5, DETECT_THRESH = 5) and specially designed to 
properly deblend the cluster galaxies (DEBLEND_NTHRESH 
= 32, DEBLEND_MINCONT = 0.005). 

Cluster members were identified by their colors, and 
their light was carefully modeled and subtracted from the 
images (Zekser et al., in prep.). This aids greatly in the 
detection of background objects, especially demagnified 
multiple images located behind the main cluster galax- 
ies. These central images provide valuable constraints 
to our lensing mass model. This galaxy subtraction 
also improves our photometry measurements and thus 
photometric redshifts for any background galaxy images 
strongly bathed in the light of a nearby cluster galaxy. 

We are now able to re-run SExtractor and detect many 
galaxies revealed by the subtraction of the cluster galax- 
ies. But the detection and object segmentation (the art 
of assigning each pixel to a given object) are still not per- 
fect. We inspect all of the object segmentations and edit 
them "by hand" where necessary. We also add a few mul- 
tiple images which are predicted by our mass model and 
are visible in the images but managed to escape detec- 
tion, either due to a b right neighbor or otherwise. Using 
the SExSeg package (Coe et al. 20061, we are able to 
force our revised object definitions into SExtractor for 
photometric analysis. 

2.3. Photometry 

With observations obtained from a wide range of tele- 
scopes, both from space and from the ground, care 
must be taken to obtain robust PSF-corrected aperture- 
matched photometry. This proves especially crucial for 
the faint lensed background galaxies we are most inter- 
ested in. Thus we use the software package ColorPro, 
which we developed and applied previously to obtain 
robus t photometry of galaxies in the UDF (Coe et al. 
20061. As in our UDF analysis, we also recalibrate the 
photometric zeropoints of several of our images. This 
procedure is described below. 



2.4. Photometric Redshifts 



LensPerfect Abell 1689 Analysis 



TABLE 3 
Abell 1689 Observations 



Telescope/Camera/Detector 



Filter 



FWHM Area Magnitude Limif 



ACS/WFC 

ACS/WFC 

ACS/WFC 

ACS/WFC 

DuPont Telescope ® Las Campanas 

Nordic Optical Telescope (NOT) @ La Palma 

Keck II / LRIS 

Keck II / LRIS 

Keck II / LRIS 

ESO NTT / SOFI 

ESO NTT / SOFI 

ESO NTT / SOFI 



g' 


O'.'ll 


11.8 D' 


27.24 


r' 


O'.'IO 


11.8 D' 


27.00 


i' 


O'.'IO 


11.8 D' 


26.92 


z' 


O'.'ll 


11.8 D' 


26.50 


U Johnson 


0'.'93 


72.3 D' 




B NOT 


I'.'IS 


41.0 D' 




VLRIS 


O'.'Sl 


80.1 D' 


27.17 


iJLRIS 


0'.'74 


93.1 D' 


26.96 


/LRIS 


0'.'68 


88.7 D' 


26.41 


J SOFI 


0'.'87 


50.4 D' 


25.12 


H Johnson 


0'.'99 


44.9 D' 


24.22 


Ks SOFI 


0'.'82 


50.4 D' 


24.06 



AB, 5-a" within a 0.8 ' diameter aperture (except for ACS, which are 10-cr, 0.2" diameter aperture) 



Photometric redshifts were obtained for the objects 
in our Abell 1689 catalog using an updated version 
1.99.2 of the Bayes ian Photometric Redshift software 
BPZ (Beni'tez 2000). This version features the recal- 
ibrated (JWVV+aB SE P (spectral energy distribution) 
templates introduced in |Benitez et al. ( |2004[) plus the tw o 
younger starburst template s added In [(Jo eetal (2006). 

In our UDF analysis (Coe et al. 200S), we reported a 
photometric redshift accuracy ot Az = 0.04(1 + Zspec)- 
Here we expect to achieve similar, or perhaps slightly 
worse, precision. The A1689 exposures are not as deep 
as the UDF, although in regions of high magnification 
(/z > 6; Am > 2) we can actually detect fainter galaxies 
than possible in the UDF. The images are very crowded 
with cluster galaxies whose light may contaminate that of 
neighbors (despite our best efforts to model and subtract 
this light) affecting the photometry and thus photometric 
redshifts. 

Abell 1689 has been the target of several spectroscopic 
redshift campai gns, includi ng Tcague et al. (1990); Fort 
(|1997); Balogh et al. (2002^; |Duc~et~ar ~ 

taT: ■ ' 



et al. (11997) 



b'Tye et 



Balogh et al. 

2007). Of the spectroscopic redshifts 



, , , , (2002, , ^ 

published in these works, 113 lie within the ACS FOV 
These are presented in Table |4] Of these objects, three 
are stars, another (#172) is half outside the ACS FOV, 
and for another (#194) it was unclear which object was 
being referenced by the published coordinates. What 
remain are 108 galaxies. In our first attempt to obtain 
photometric redshifts for these galaxies, it was clear that 
our SED templates provided a poor fit for the observed 
U and B magnitudes. 

To reveal flux miscalibrations in any of the filters, we 
apply techniques sim ilar to that used in our study of the 
UDF ( Coe et al.|200 6). We select cluster ellipticals based 
on their spectroscopic redshift, BPZ spectral type fit, and 
visual confirmation in the ACS images. We then use BPZ 
to fit SEDs to the photometry of these objects, fixing 
the redshifts at the spectroscopically-determined values. 
We find the following offsets between the observed and 
predicted magnitudes: V -0.08, R -0.10, / -0.00, J -0.08, 
H -0.06, Ks -hO.Ol, g' -0.15, r' -h0.03, i' -0.04, z' -^0.02. 
We subtract these biases from our measured magnitudes. 
Our U and B-band images did not arrive well calibrated, 
so we used this procedure to calibrate them as well. 

Given our recalibrated fiuxes, we rerun BPZ on all our 
galaxies without constraining the redshifts to the correct 
values as above. We measure of goodness of fit Xmod 
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Fig. 1. — Example of an excellent BPZ SED fit to the photom- 
etry of a cluster elliptical (#198), after flux recalibrations of all 
filters (see text). The photometric redshift is correctly derived as 
z = 0.18 it 0.12 (95% confidence). The colored circles give the ob- 
served AB magnitudes with uncertainties in 12 filters, while the 
blue rectangles give the SED model predicted magnitudes (with 
rough uncertainties) in those filters. The colors serve to distin- 
guish the telescopes / filters: purple - Las Campanas, La Palma 
(UB); green - ACS {g'r'i'z'); orange - Keck/LRIS (VRI); red - 
ESO NTT/SOFI (JHKs). 
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Spectroscopic Redshift 

Fig. 2. — Bayesian photometric redshifts plotted versus 108 spec- 
troscopic redshifts obtained within the ACS FOV (TableH. Good 
BPZ SED fits (Xmod < ^) ^'^^ plotted in blue, while poorer fits 
to the photometry are plotted in yellow. For the 89 galaxies with 
good SED fits, the photometric redshifts agree with the spectro- 
scopic values to within A^ = 0.07(1 -I- Zgpgc)- If we recursively 
remove outli ers with m ore than three times the RMS deviation Az 
(as iU jCoe et al.|2006| , two galaxies are pruned and the rest agree 
to witmn Z\2 = U.l 
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between observed and model fluxes with some uncer - 
tainty assigned to the model fluxes (see Coe et al.|2006 ) . 
An example of an excellent {Xmod ~ 0.03) and correct 
(Az = 0.01) SED fit to a cluster elliptical obtained with 
recalibrated fluxes is shown in Fig. fl] Of the 108 galax- 
ies with spectroscopic redshifts in the ACS FOV, we ob- 
tained good SED fits ixl^od < 1) for 89. For these 89 
objects, our Bayesian photometric redshifts agree with 
the spectroscopic redshifts to within Az = 0.07(1 + Zgpec) 
(see Fig. pi) . And if wc recursively remove o utliers with 
more than three times this deviation (as in Coe et al. 
2006 1 , then two galaxies are pruned and the remaining 



ous systems with thick multiple images as in, say, CL0024 
(jTyson et al.||1998[ [Zitrin et al. ||2009b[). Only after care- 



87 agree to within Az = 0.06(1 + Zspec)- The two outliers 
are hardly catastrophic, having Az = 0.30(1 -t-z^pec) and 
Az = 0.21(1 + Z.peJ- 

Finally, we obtain photometry and photometric red- 
shifts for our multiple images. The results are given in 
Table [5] For each redshift, a 95% confidence interval 
is quoted as well as Xmod- Note that in B05 the BPZ 
prior was modified. The redshift range z < 0.7 was ex- 
cluded and the prior allowed for observed fluxes up to 
20 X magnified. We have not implemented this modified 
prior here. 

Note that these redshifts were obtained using photome- 
try from all 12 filters. We have also obtained photometric 
redshifts based on photometry obtained in the four ACS 
images with the light of cluster galaxies modeled and sub- 
tracted. We are unable to apply this galaxy subtraction 
in a consistent and robust manner to the ground-based 
images due to their significantly worse PSFs. Yet even 
with four filters, we obtain improved photometric red- 
shifts for those objects whose light is significantly con- 
taminated by nearby cluster galaxies. 

3. MULTIPLE IMAGES 

Building on previous work, we present a catalog of 135 
images of 42 background galaxies. We contribute 20 new 
candidate multiple images of 8 galaxies, along with a 
few tweaks to previous identifications. We discard three 
central demagnified images in use since B05, finding these 
identifications suspect; our technique is more sensitive to 
the p ositions of central images than other techniques (see 
^3.2). We add one new central image identification. 

We also identify multiple knots in 42 of the galaxy im- 
ages which we use as additional constraints. Constrain- 
ing the positions of three non-coUinear knots in an im- 
age is equivalent to constraining the precisely-measured 
shear and magnification of that image. Our mass map so- 
lutions perfectly reproduce the observed positions of 168 
multiple images of 55 knots residing within 135 images 
of 42 galaxies. The final set of multiple images used in 
this work is given in Table [5] and shown in Fig.[3j Close- 
ups of all the images are shown in Fig. |16[ Aaditional 
knots are visible in these images and listedm Table |6] In 
this section we discuss the multiple image identification 
as well as the redshifts for these systems. 

3.1. Multiple Image Identification 

The original analysis of the ACS A1689 images (B05) 
yielded 106 multiple images of 30 background galaxies. 
This was a truly pioneering effort as the first multiple im- 
age identifications were the most difficult. The relatively 
steep mass profile of A1689 near the Einstein radius pro- 
duces relatively thin arcs. Thus there are no truly obvi- 



fui study of the color image did B05 discover the image 
system 1-2, a pair of pale green and blue specks that re- 
peats five times about the image, leaving no doubt as 
to its identification. Once these first multiple images are 
identified, an initial mass model may be obtained, greatly 
facilitating the identification of further image systems. 

B05 identified many image systems, but they did not 
attempt to identify all. Additional systems have since 
been proposed, and we propose still more in this paper. 
Of the 30 image systems proposed by B05, doubts have 
been since been raised about three. We further call into 
question three of their central image identifications. 

Z06 (their §8.4) experimented with excluding some im- 
age systems from their modeling on the grounds that 
they yielded larger than average scatter in the delensed 
positions. But they stopped short of calling any of B05's 
systems into question, as large scatter may also result 
from shortcomings of the mass model. Among the prob- 
lem systems they cited were systems 1-2 and 15 which 
have ironclad identifications as multiple images from vi- 
sual inspection. Images near critical curves such as 
these often prove difficult to fit for conventional modeling 
methods, as the predicted image positions are unstable 
with respect to small variations in the mass model. For 
this reason, L07 also found it necessary to remove the 
giant arc (system 8) from their modeling. 

Some more systematic changes to the B05 catalog were 
made by H06 and L07. They discarded and altered a few 
systems which appear to have been misidentified in the 
original B05 analysis. H06 (see discussion in their §3.3) 
split B05's system 12 in two: a new streamlined system 
12 comprised of just 12b and 12c; and a new 4-image 
system 13 comprised of 12a, 12d, and two new counter 
images. (Note this offsets their numbering relative to 
B05 for B05's systems 13-19.) Splitting 12a and 12b 
into different systems was a bold proposition, given that 
they were known to have the same spectroscopic redshift 
of z = 1.83. But L07 concurred with this assessment 
(see their §4.2), and so do we. We adopt L07's conven- 
tion of assigning the number 31 to H06's system 13, thus 
preserving the B05 numbering system. 

B05 may have misidentified the counterimages of the 
fold arc 12bc, but this arc almost certainly has other 
counterimages. We identified three counter images 12fgh 
plus an image 12i which is a counter image of 12g due to 
strong galaxy-galaxy lensing. (Note these designations, 
12fghi, were intended to avoid confusion with B05's orig- 
inal identifications of images 12abe. We did not intend to 
suggest this is a 9-image system. There are but 6 images 
in our new system 12.) 

H06 also discarded B05's 3-image system 20 located in 
the top-left corner of the image at large radius. Not only 
do the images fit poorly with the mass model, but the im- 
ages show slightly different morphologies. L07 concurred, 
adding that the spectra are also somewhat different. We 
did not attempt to add this system to our mass model. 

L07 also took issue with B05's system 25, finding a dif- 
ferent counterimage to 25b and naming the new system 
33. We concur, finding the new system a much better 
fit. As for B05's 25a (a greenish radial arc), we identify 
a new counterimage candidate at fairly large radius, and 
we rename the system 45. 
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Fig. 3. — Multiple images in Abell 1689 labeled on the STScI 3.2' X 3.3' ACS g'r'z' color image. Most of the labels are directly above 
the images, but some have been offset slightly to minimize clutter. The exact coordinates of each image are given in Table [s] A compass 
provides the orientation (a 115° roll angle), while the lengths of the arrows provide the scale (20" ~ 62 kpc). 



H06 identified one brand-new system (their 31), a small 
pair of blue images along the fold inside the main sub- 
clump. L07 "rediscovered" this system, naming it 36. 
L07 also propose the following brand-new systems with 
numbers of images given in parentheses: 32 (4), 35 (3), 
36 (2), 40 (2). We include all of these in our mass model. 

We have proceeded to identify 8 new candidate sys- 
tems containing a total of 20 multiple images. We have 
named these (with number of images given): 41 (3), 42 
(4), 44 (2), 45 (2), 46 (2), 48 (2), 49 (2), 50 (3). As men- 
tioned above, we recycled B05's image 25a for use in our 
system 45, but the other 19 multiple images are brand 
new identifications. 

We do not have spectroscopic redshifts for our new 
systems; thus our two-image systems do not currently 
provide strong constraints. In the course of our mass 
modeling, we add our new systems at the end and find 



that our two-image systems have little effect on the mass 
map. 

We also mention in passing H06's "system" 32, a sin- 
gle large but low surface brightness arc just outside the 
image pair 6cd in the main subclump. They were unable 
to identify a counterimage for the arc and neither were 
we. We do not utilize this arc in our modeling. 

3.2. Demagnified Central Images 

The identification of central images is crucial to con- 
straining the inner mass profile of any lens. This is gen- 
erally a difficult task both because central images are 
strongly demagnified and because this faint light is over- 
whelmed by the bright galaxy or galaxies which make 
up the lens. We have carefully modeled and subtracted 
most of the cluster galaxies from each of the g' , r', i' , 
and z' ACS images (Zekser et al., in prep.). The re- 



Coe et al. 



combined galaxy-subtracted color image reveals many 
colored specks in the central region of the cluster. A 
great number of these specks a re globular clusters be- 
longing to A1689 (Mie ske et al. 2004, 2005). The rest 
may be demagnified multiple images. Asidelrom sorting 
out this confusion, we also have to contend with resid- 
uals which persist from the galaxy subtraction. Correct 
identification of a central multiple image relies on finding 
the proper color speck (or group of specks for paired im- 
age systems) at approximately the location predicted by 
the model. Often we find several specks of approximately 
the correct color in approximately the predicted location. 
Thus it should not be surprising when we claim that a 
few mistakes may have been made previously. Conven- 
tional mass modeling has allowed these mistakes to go 
unnoticed by BOS and in subsequent studies. 

As mentioned above, conventional modeling methods 
may be very sensitive to the positions of images in re- 
gions of high magnification. Conversely, they generally 
will not be very sensitive to the exact positions of de- 
magnified images. Our model-predicted positions for the 
central images 7c, 8e, and 19e are "only" offset by ~ 2'.'5, 
7", and 2'.'5 from the B05 positions, respectively. Thus 
these offsets may be easily missed by routines that min- 
imize offsets in the image plane. (Note that attempts 
may be made to normalize the offsets by the local mag- 
nification, but this is often not attempted as it can lead 
to instability in the optimization routine.) If optimiza- 
tion is instead performed in the source plane (technically 
inferior but much quicker), the method might actually 
be more sensitive to the offsets of central images. L07 
do find higher than average (0'.'4) offsets in the source 
plane for systems 7 (1'.'09) and 19 (^.'48). And as men- 
tioned above, L07 found such large errors for system 8, 
that they excluded it from their analysis. The incorrect 
identification of 8e, 7" from the position we predict, may 
have contributed to their errors as much as the unsta- 
ble model positions of 8ab (near the critical curve for 
that redshift). The high offset for the 3- image system 7 
may have raised alarms in L07's analysis. But they find 
similarly large offsets for the pair of systems 26-27, and 
we find no quarrel with this pair. The 26-27 pair does, 
however, strongly require asymmetry in the center of the 
mass map, which had not been observed by previous au- 
thors. In conventional mass modeling, it is difficult to 
say when a larger than average offset is a misidentifica- 
tion and when it may simply indicate a shortcoming of 
the mass model. 

LensPerfect, on the other hand, is extremely sensitive 
to the relative positions of all of the central images. For 
example, a central image incorrectly identified to the 
wrong side of another central image is generally disas- 
trous for the mass model, causing the deflection field to 
get tangled in itself. 

We reexamined all of BOS's central image identifica- 
tions. The pair lf-2e is confidently identified by the im- 
ages' colors and proximity to one another. As for the 
rest (4e, 7c, 8e, 19e, 10-15-18c, lie, 22-23b, 26-27c, 32d, 
35cPM we purged them all from the image list and ob- 
tained a mass model solution without them. We then 

^^ Dashes link objects (e.g., 10, 15, and 18) assumed to be in 
physical groups with redshifts constrained to be equal in our model. 
We are not referring to objects 10 through 18, inclusive. 



re-added the multiple images one by one to our model. 
We found that a few central images (7c, 8e, 19e) did 
not fit well with the rest, producing aphysical models 
when added (even after all of the source positions were 
allowed to shuffle to new positions in search of a physical 
model). We are unable to securely identify replacement 
central images for systems 7, 8, and 19, as there are too 
many similarly colored specks in the area which confuse 
the issue. 

Red demagnifled central images stand out more and 
are easier to identify. We identify one new central image 
candidate 28c which we incorporate into our mass model. 
It flts easily into our model, not signiflcantly affecting it. 

3.3. Additional Knots 

With LensPerfect we obtain mass map solutions which 
delens the input centroids of all images of a given system 
to the exact same position in the source plane. But this 
alignment does not guarantee that the delensed images 
will have the same shape or orientation. Thus we identify 
additional knots where possible in the multiple images. 
These knots are labeled in Fig. [16] If the delensed po- 
sitions of these knots do not align well naturally given 
our mass model, then we add them as constraints and 
force them to align. In Table [6] we give the positions of 
these additional knots which we have constrained in our 
mass model (in addition to the centroids / primary knots 
listed in Table Isl) . These knots are also labeled in green 
in Fig. [T7| In all, our model incorporates 168 observed 
(lensed) positions of 55 knots. 

These additional constraints further improve the ac- 
curacy of our mass map (assuming our multiple image 
identiflcations are robust, as we believe they are). Con- 
straining three non-collinear knots in an image effectively 
uses both the obs erved shear and re lative magniflcation 
of that image. In Coe et al. (2008) we compared mass 



maps of MS1358 given a single multiply-imaged galaxy 
both with and without additional knots constrained. The 
additional knots add signiflcant detail to the mass map. 

3.4. Redshifts 

In the original B05 analysis, spectroscopic redshifts 
were available for 5 systems: lad 3.04, 2a 2.54, 7a 4.87, 
10a 1.37, and 12ab 1.82. Since then, system 12 has been 
split in 2 by II06 (our 12 and 31) and the redshift of 10a 
has been called into question by L07, with the net result 
being that we still had 5 systems with available redshifts. 
Our initial analyses made use of these spectroscopic red- 
shifts alone. 

L07 have since contributed spectroscopic redshifts for 
another 19 systems (4a, 5ac, Bab, 10-15-18a, 11a, 14a, 
17c, 19d, 22a, 24c, 29d, 30c, 32b, 33a, 35a, 36a, 40a), 
bringing the total to 24 systems with spectroscopic red- 
shifts (Table [5]) . We have compared our results before 
and after incorporating these redshifts into our mass 
modeling. We find the substructure shifts somewhat, 
but overall the mass models appear to be very similar 
qualitatively. Our "optimized" redshifts (described be- 
low) were generally close to the spectroscopic redshifts, 
with a scatter of Az « 0.06(1 -I- Zg) after pruning one 
outlier. 

For those systems without spectroscopic redshifts, we 
use photometric redshifts (Q as initial guesses in our 



LensPerfect Abell 1689 Analysis 



optimization routine (q4|. We allow these redshifts to 
wander but they incur a penalty for doing so. A rough 
uncertainty {Az = 0.5,0.7,1.0,2.0) is assigned to each 
redshift according to our relative confidence in it, and 
deviations from the input values are divided by these un- 
certainties. Given these normalized deviations, we then 
take the RMS and add this to our penalty evaluation. As 
we are modeling 43 image systems, the total RMS will 
not be very sensitive to single outliers. This approach 
will allow individual redshifts to wander unacceptably 
far. In order to "leave no redshift behind," we also find 
the maximum normalized deviation and add this to our 
penalty function. Our input redshifts and uncertainties 
along with final optimized redshifts for each system are 
given in Table [7J 

We could attempt to penalize redshift deviations more 
scientifically by using the redshift probability distribu- 
tions P{z) returned from BPZ for each object. However, 
photo-z uncertainties are ofte n underestimated by cur - 
rent methods including BPZ (Hildebrandt et al. 2008). 
Thus we prefer to assign simple and rather generous un- 
certainties to the redshifts and allow them to naturally 
obtain their optimal values based on the mass model. 

Some sets of images appear to be physically linked: 
10-15-18, 22-23, 24-29, 26-27. During our optimization 
procedure, we find that within each of these sets, all 
of the redshifts gravitate toward common values (10-15- 
18: 2.00-2.14-1.96; 26-27: 1.98; 24-29: 1.91, 22-23: 1.4, 
1.46). Thus we take the liberty of fixing all of the red- 
shifts to be equal within each set. For example, when the 
redshift of object 10 is optimized, the redshifts of objects 
15 and 18 are forced to follow. We believe these systems 
to be physically linked but any one of them might instead 
be a chance alignment. Object 15, for example, did ob- 
tain a slightly higher redshift than 10 and 18 when left 
free as noted above. 

Objects 1-2 also appear to be physically linked, but 
our models say otherwise. The objects do have different 
spectroscopic redshifts, but that obtained for object 2 
is somewhat uncertain, as its spectra appears to show 
absorption from three separate system s along the line o f 
sight at redshifts 2.53, 2.87, and 3.04 ( |Frye et al.| [20Cl7t. 
We tried fixing the redshift of object 1 to 3.04 and giving 
object 2 an initial guess of the same redshift {z = 3.04) 
but allowing it to wander (incurring a small penalty for 
doing so). A lower redshift of ^ 2.5 was clearly preferred 
for system 2, thus s upporting the spe ctroscopic redshift 
of 2.53 preferred by|Fryeet aL|(|2007[). 



3.5. Coordinate System 

The final catalog of multiple images used in our mod- 
eling is given in Table [5) Wc provide coordinates in both 
(RA, Dec) and in (x, y). Our (x, y) coordinate syst em is 



based on the original APSIS ( [Blakeslee et aL|2003[ ) ACS 
GTO pipeline reductions. These images are 4421 x 4525 
pixels (1 pixel = 0.05"). Our bottom-left pixel is cen- 
tered at (1,1) as in SExtractor and ds9. North and East 
are 115° clockwise from up and left, respectively. 

Based on the central 3853 x 4000 area of these images, 
STScI released a g'r' z' color imageFj An offset of (350, 
232) may be subtracted from our coordinates to obtain 
coordinates in the color image. The ACS images were 



later reprocessed by APSIS yielding images 4379 x 4481 
pixels. The improvements included better correction for 
the geometrical distortion of the images. Nevertheless, 
the offset in coordinates between the original and sub- 
sequent processed images is a nearly constant (21, 22) 
across the entire image, only deviating by a pixel in x in 
the top-left and bottom-right corners. 

4. LENSPERFECT MASS MAP RECONSTRUCTION 

LensPerfect is a novel approach to gravitational lens 
mass map reconstruction. The 100+ strong lensing fea- 
tures produced by Abell 1689 present us with a large 
puzzle. We must produce a mass model of A1689 with 
the correct amounts of mass in all the right places to 
deflect light from 30-1- background galaxies into multi- 
ple paths such that they arrive at the 100-f positions 
observed. 

Most strong lensing analysis methods construct many 
possible models and then iterate to find that which best 
matches the data. LensPerfect instead uses direct ma- 
trix inversion to find perfect solutions to the input data. 
Using LensPerfect, we may, for the first time, obtain a 
mass map solution which perfectljP^I reproduces the in- 
put positions of all 100+ multiple images observed in 
Abell 1689. 

LensPerfect makes no assumptions about light trac- 
ing mass. Non-LTM models a re common in analyses of 
WL or coinbined SL-F-WL (e.g., Bradac et al. 2006; Diego] 



et al.|2007[|Deb et al. 2009 ). Dedicated SL analysis meth^ 
ods are able to process greater numbers of multiple im- 
ages. Non-LTM SL analysis metho ds inclu de PixeLens 
jSaha fc Winiams|2004|[Coles|2008l ), SLAP (pegoeTaL 



2005a b), and m ethods developed by Liesenborgs et aT? 
t2006[ 12009] ) and [Julio fc Kneib] ( |2009[ ). ('I'he latter in- 



cludes both LTM and non-L'i'M components.) We note 
non-LTM methods are also used in SL modellin g of ex- 



tende d images lensed by individual galaxies (e.g., Vegetti 



et al.|2 009 

LensPerfect was made possible by a recent advance 
in the field of mathematics (Fusclier 2006 2007) . The 
method was described in detail in jCoe et al.| ( |200g| ). Here 
we provide a brief outline of the procedure. 

Image defiection by a gravitatio nal lens is governe d 



by a few simple equations (e.g., Wambsganss 1998) 
Gi ven the bend an gle a of light due to mass predicted 
by Einstein (1916), we can derive the defiection of light 
due to a mass sheet with surface density ac as a function 
of position 6 in the lens / image plane: 



d{0) = - / (fe'nie') 



with the simple corresponding inverse relation: 



V -a = 2k. 



(1) 



(2) 



http://hubblesite.org/newscenter/archive/releases/2003/01/ 



The surface density k = Y,/Yicrit is defined in units of 
the critical density at the epoch of the lens. The critical 

^^ Again, while the image positions do have (small) observational 
uncertainties of a pixel (0.05") or so, our solutions do perfectly fit 
the data as input. We may vary the image positions within this 
uncertainty to produce other valid solutions, however this is not a 
significant uncertainty in our models. 
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density is that generally required for multiple images to 
be produced. It is a function of source redshift as given 
by: 



^crit — 



Ds 



AttGDlDls' 



(3) 



involving a ratio of the angular-diameter distances from 
observer to source Dg = Dji{0,zg), observer to lens 
Dl = Da{0,zl), and lens to source Dls = Da{zl,zs)- 
For a flat universe (f2 = Qm + ^A — 1), angula r- diamet er 
distances are calculated as follows {Fukugita et al.|[l992 
filled beam approximation; see also iHoggiil999 1 : 



Da{zi,Z2) = 



dz' 



I + Z2J,, H{z'y 



where the Hubble parameter varies with redshift as: 



(4) 



(5) 



(Formulae for non- flat cosmologies can be found in Coe 
fc Moustakas||2009 for example.) 

l^'rom simple geometry we find that the deflection 
angle a (and thus the critical density Ticrit) is a function 
of redshift. The deflection is greatest (aoo) for a source 
at infinite redshift. For sources at less than infinite 
redshift, this deflection is reduced by the distance ratio: 



(Dls\ . 



(6) 



Thus the problem of mass map reconstruction can be 
reduced to determining the deflection fleld with all de- 
flections scaled to a common redshift (e.g., c?oo), at which 
point we simply take the divergence and divide by 2 
to obtain the mass map (Eq. [If. The deflection field 

d{6) = — j3 may be measured at the multiple image po- 
sitions 9 once source positions j3 are determined. How- 
ever, in order to take its divergence, the deflection field 
must be solved for as a continuous function of position 
(or at least defined on a regular grid). Our in terpolated 
deflec tion fleld must also be curl-free (see e.g., Coe et al. 
2008| . 

Only recently were the mathematical tools developed 
that enable us to obtain a curl- free interpolation of a 
vector fleld (Fuselier 2006 2007). The technique uses 
direct matrix inversion to obtain a solution which ex- 
actly matches the vectors at the given data points and 
interpolates or extrapolates it elsewhere. The solution is 
composed of radial basis function (RBFs) each placed at 
the position of an observed multiple image. Each basis 
function has two free parameters (amplitude and rota- 
tion) equal to the number of constraints (the x and y 
coordinates of the image) . While some non-LTM meth- 
ods have many more free parameters than constraints, 
ours does not. 

In gravitational lensing, our vector (deflection) fleld is 

not defined until we assume source positions /3 for our 
lensed galaxies. As we add each galaxy to our model, 
we can obtain a good initial guess for each /3. We then 
perturb all of our source positions. Each arrangement 



of source positions yields a new mass map. Iterating 
over various arrangements of source positions, we find a 
range of mass maps all of which perfectly reproduce the 
observed image positions. Among these, we select the 
"most physical" mass map using a set of non-restrictive 
criteria. Aside from the requirement that the mass map 
be po sitive, these criteria (described in detail in Coe et al. 
2008 1 promote mass maps which are smooth, decrease 



outward from the center on average, and are azimuthally 
symmetric (small scatter in radial bins, with extra penal- 
ties for "tunnels"). 

Lensing generally constrains the projected mass within 
the Einstein radius, or more precisely, the region within 
the multiple image positions. This region is known as 
the "convex hull" in the language of our interpolation 
scheme. Outside the convex hull, our solutions are ill 
defined and in fact drop off to zero (and even negative 
values) too quickly. Our mass models should generally 
be disregarded outside this region 



Since pub lishing the LensPerfect method paper (Coe 
et al. 2008), we have made small changes in the ex- 
act implementation of these penalties and their relative 
weights. And where (for the purposes of calculating 
penalties only) we had evaluated the mass map on a 
41 X 41 grid we now evaluate it on a finer 81 x 81 grid 
within the convex hull and a coarser 21 x 21 grid outside. 
(With the 41 x 41 grid, we found one particular "tunnel" 
was escaping detection.) Finally, we have adde d th e red- 
shift penalty function as described already in §3.4| 

We emphasize that our mass models are not "grid- 
based" . The radial basis functions are instead placed at 
the positions of the multiple images, as described above. 
The resulting mass model has a smooth functional form 
and can be calculated at any desired coordinates. How- 
ever in order to evaluate and present the results, we gen- 
erally calculate the mass model on a regular grid. 

5. MASS MODELS 

Here we present non-LTM mass models which perfectly 
reproduce the observed positions of 168 multiple images 
of 55 knots within 135 images of 42 galaxies, strongly 
lensed by A1689. We stress that there is no unique solu- 
tion, and we do obtain a range of solutions which allow 
us to estimate our uncertainties (j |5.2[ ). However first we 
present the most "physical" solution found by our op- 
timization scheme described above. This optimization 
took two weeks to run on a MacBook Pro laptop. The 
process runs quickly at first (a few minutes per galaxy 
added) but slows as more galaxies are added. Galaxies 
were added in the order presented in Table [7j 

5.1. Most Physical Mass Model 

In Fig. [4] we present our most physical strong lens- 
ing mass model of A1689. The mass map contours are 
laid over the ACS STScI g'r'z' color image in Fig. [5| 
Our model is constrained best near the multiple im- 
ages (shown in pink), interpolated between them, and 
is highly uncertain where it is extrapolated outside the 
outermost multiple images (our "convex hull"), traced in 
white in Fig. |4] and black in Fig. [5] In fact our model 
falls off to zero too quickly outside the convex hull. 

From inspection of Fig. [5] we find that our mass model 
resolves halos of perhaps ten or so galaxies which are 
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Fig. 4. — Mass map solution for Abell 1689 which perfectly re- 
produces the 168 observed (strongly lensed) positions of 55 galaxy 
knots / centroids. Mass is plotted in units of Koo (critical den- 
sity for a source at Zg = oo) and cut to the same 3.2' X 3.3' field 
of view as the STScI color image (previous and following figures). 
The black line indicates the convex hull. Multiple images are found 
within this region. Outside this region, our solution is highly un- 
certain and in fact falls off to zero too quickly. A red line marks 
the K = contour. We stress that this solution is not unique, but 
had the highest "physicality" of all solutions we explored. Angular 
diameter distances are given along the axes. 

members of A1689. Although we have made no assump- 
tions about LTM, much of our mass model's substruc- 
ture does coincide with luminous galaxies. Determining 
masses for these individual galaxy halos would be useful 
though far from trivial, as we discuss in ^\ 

There are, however, some potentially interesting off- 
sets between mass clumps and luminous galaxies. For 
example, the mass near the BCG appears more pinched 
than the distribution of galaxies just above and to the 
left. Several multiple images (plotted as pink squares) 
are present in this region lending confidence to our mass 
model there. Furthermore, when we tested our method 
on simulated len sing with a similar mass distribution 
(Coe et al. 2008), our recovered mass maps exhibited 



no such pinching 

More quantitative and robust conclusions about the 
alignment of mass and light will await future work (§7]). 
We will perform tests to determine how robustly sub- 
structure clumps are identified and their positions deter- 
mined. 

We estimate that our mass map resolves substructures 
~ 23 kpc across within Re ~ 150 kpc of the core (an- 
gular diameter distances). This estimate is based on the 
density of multiple images, N(P = 7ri?|;, where N = 135 
images, and we find the average separation among these 
to be d « 23 kpc. Each multiple image provides a con- 
straint on our deflection field and thus our mass model. 
The observed multiple images are more densely packed 
in some regions; thus the mass map resolution will be 
greater there and lesser elsewhere. 

This is the highest resolution mass map to date of any 
galaxy cluster without assuming LTM. [Leonard et al.| 
(2007) pres ent a SL-I-WL m ass n iap with pixels ^ 65 



kpc across. |Saha et al.| ( |2006[ ) and [Julio fc Kneib| ( |2009 1 
use smaller pixels but computational issues limit theni 



to ~ 30 multiple images per solution. We estimate their 
effective resolution to be ~ 50 kpc, or about one-fourth 
our 2-D resolution (half along each axis). By using four 
times the number of constraints (multiple images), we 
obtain a mass model with four times the resolution. 

Finally in Fig [17] we show our multiple images as de- 
lensed to the source plane by our mass model. Note that 
each constrained knot in each system is delensed to the 
same source position in each multiple image. To be clear, 
we do properly model extended images as the multiple 
knots in each image map back to different locations in the 
source plane. For example, in the large arc (Sab), three 
distinct knots are identified and constrained in each im- 
age (a and b) . Knots 8a0 and 8b0 map back to one point 
in the source plane, 8al and 8bl map back to another, 
and 8a2 and 8b2 map back to a third point in the source 
plane. 

5.2. Mass Model Ensemble and Uncertainties 

The mass model solution presented above is not 
unique. Using the optimization procedure described in 

fwe explored a wide range of source positions and red- 
fts, and found a set which produced this most "phys- 
ical" mass map according to our criteria and without 
letting the redshifts stray too far from their input val- 
ues. 

We estimate the uncertainties in our mass model by 
exploring an ensemble of mass model solutions. In future 
work we will develop algorithms to explore this solution 
space more thoroughly, as we describe in 37] 

Here we settle for a proxy ensemble of mass models, 
based on a broad but non-exhaustive search of our pa- 
rameter space. This ensemble consists of 54 solutions 
which we obtain as we add galaxies building up to our 
final "best" solution. The first model in our ensemble 
is optimized given systems 1 and 2 only, and our last 
includes all 55 knots. 

It is unclear whether this technique should be expected 
to overestimate or underestimate somewhat our actual 
uncertainties. We include too broad a range of solutions 
by including those that only fit some of the data, yet we 
may not be thoroughly exploring the solution space. 

This technique does capture some of the systematic 
uncertainties which would result from adopting various 
subsets of the multiple image systems. Previous authors 
have done just that: used subsets of our multiple im- 
ages. In previous work there has also been some varia- 
tion in these identifications. One or more of the systems 
presented in our work may yet prove incorrect, and we 
capture some of those uncertainties here. 

Our estimated uncertainties on the radial mass profile 
appear to have the correct form. Mass enclosed within 
the Einstein radius Re ^ 47" i s con strained more tightly 
than mass within other radii {{ 6.1 ). 

The exact structure we resolve in the center is sensi- 
tive to the identifications of demagnified central images, 
which can be fairly uncertain. However we believe our 
modeling method enables robu st identification of central 
images. As discussed in §3.2[ we have reexamined the 
B05 central image identifications, purged those which do 
not fit well with the rest, and identified a new central 
image candidate (28c) which we include in our model. 

In Fig. [6] we plot the critical curves for sources at 
redshifts Zs = 1,2, and 7. The exact shapes of the critical 
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Fig. 5. — Mass map contours in units of Kqo = 1/3 laid over the 3.2' X 3.3' STScI ACS g'r' z' color image. The outermost contour, 
Koo = 0, was also plotted in the previous figure. Pink squares indicate the 135 multiple image positions all perfectly reproduced by our 
model, and the white line indicates the convex hull. Outside this region, our solution should be disregarded. This solution is not unique 
but was the "most physical" we found. 



curves are not well constrained by the data alone. They 
are sensitive to the exact distribution of substructure in 
the mass map, which for this purpose is not resolved 
sufficiently by our models given the number of multiple 
images. The critical curves would be better constrained 
by additional multiple images or a well-informed prior on 
the mass distribution. Even in this limiting case of broad 
model freedom, we find that points along our Zg — 7 
critical curve can expect to have // > 6 68% of the time. 
We remark that in general, non-LTM methods proba- 
bly overestimate uncertainties while LTM methods prob- 
ably underestimate them somewhat. LTM methods do 
not explore the full range of solutions (including asymme- 
tries and other deviations from LTM) which may repro- 
duce the data. Furthermore, they are sometimes forced 
to discriminate between a "best" solution which repro- 
duces the data at, say, 50-cr (2.5"), and "deviant" solu- 



tions which reproduce the data at, say, 60-cr (3.0"). Non- 
LTM methods, on the other hand, may include, at worst, 
a broad range of unphysical solutions or, at least, solu- 
tions which do not take advantage of the strong obser- 
vational priors available, namely the observed positions 
of the lensing galaxies. An ideal method would use LTM 
as a prior while allowing for deviations (aTl). This prior 
might be referred to as "LATM" , or light approximately 
traces mass. 

6. MASS PROFILE AND CONCENTRATION 

As discussed in ^ recent studies show that galaxy 
clusters formed in nature may be more centrally concen- 
trated than their counterparts formed in simulations. In 
this section, we obtain new estimates for the mass con- 
centration of A1689. Our results support previous anal- 
yses which claim Al689's concentration is higher than 
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Fig. 6. — Left: Critical curves for sources at redshifts Zs = 1, 2, and 7 laid over the 3.2' X 3.3' STScI ACS g'r'z' color image. The 
loops in the curves are probably too broad, an artifact of the insufficiently resolved substructure clumps in our models. (Compare the 
curves with the pink squares indicating the multiple image positions which provide model constraints.) Given additional multiple images 
(or a well-informed prior such as LATM), finer substructure would be resolved, yielding tighter and more precise critical curves. Right: 
The great freedom in our models allows for broad variation in the critical curve shapes. We plot Zg = 7 critical curves for all solutions in 
our ensemble, from that which fits two multiple image systems (dark red) to our final solution (shown in white) which fits all 55 systems, 
including multiple knots per galaxy. The variation in magnification is such that points on our final critical curve should expect to have 
fi > 6 68% of the time. This would increase for a better constrained model. 
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Fig. 7. — Our mass map fit to NFW and Sersic profiles. For all 
points within our convex hull, we plot in light blue k.{zs = oo), 
the projected surface mass density in units of critical density for a 
source at infinite redshift. Binned averages are plotted with uncer- 
tainties measured as the scatter among an ensemble of solutions. 
Our best NFW fit is plotted as the red line. The mass profile 
obtained by BOS is also plotted for comparison. 

expected. We derive these concentration estimates from 
fits of NFW profiles simultaneously to our SL mass model 
and WL data published elsewhere. These profiles do si- 
multaneously fit both the SL and WL data well, however 
there is a slight tens ion between the two as described in 
§6.2[ First, in §6.11 we examine our SL mass profile in 
detail. 

6.1. Profile fits to the strong tensing (SL) data 
Clusters formed in simulations have mass profiles gen- 



( [Coe|[2010 ). 
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2010 1 . I'hese profile fitting functions are 



described in Appendix ffl with further details given in 
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Fig. 8. — A1689 projected mass within a cylinder as a function 
of radius R. The mass profile of our best model is plotted as the 
dashed line while the shaded region gives the mean and RMS of 
the ensemble. Note the mass is especially well constrained within 
R ~ Re- The radii of the multiple images are plotted along the 
bottom. Those along the convex hull are plotted in red. Four 
central images are outside the range of the plot at small radius. 

In Fig. [7] we plot our projected mass density radial pro- 
file K,(K) with our origin defined as the location of our 
density peak which corresponds to the location of the 
BCG. We plot the mass profile from B05 for comparison 
and find good agreement even though the analysis meth- 
ods are very different (including LTM vs. non-LTM). For 
a similar recent comparison, see ( Zitrin et al.||2010b ). 

Overplotted in Fig. [t] are our best iNl-'W (r^ =338 
kpc, C200 — 7.6 \cyir = 9.6]) and Sersic (i?e = 692.8 
kpc, Ke = 0.1007, n = 2.148 [&„ = 3.968]) fits to 
the SL data alone. The two fits track each other very 
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closely, although the more flexible Sersic profile does al 
low for a slightly shallower central slope. We could al 
low the central slope of our NFW profile to vary with a 
three-parameter "g eneralized NFW profile" (|Zhao||1996" 
Wyithe et al.||2001), but we do not explore this here. 



Kather than fitting to k{R) as is common, we can re- 
duce the uncertainties of our NFW fit parameters by fit- 
ting to M(< R). The former is projected mass surface 
density at a given radius R, while the latter is total pro- 
jected mass enclosed within a cylinder of radius R. 

We plot M(< R) from our SL modefing in Fig.ls] We 
also plot the uncertainty in M(< R) from our ensemble 
of models. This quantity, M(< R), is what lensing con- 
strains best, especially M(< Re), the mass within the 
Einstein radius. Though our models vary, all give a very 
consistent amount of mass within Re ~ 47". 

We derive an Einstein radius of Re = 47.0" ± 1.2" 
(143^4 kpc) for a lensed source at z^ = 2. This value 
increases with redshift to Re ^ 52" for a source at Zg = 
7. No cluster is perfectly symmetric, so the "correct" 
definition of Einstein radius is a bit ambiguous. As done 
elsewhere, we find that radius within which (k) — 1, 
the average mass surface density is equal to the critical 
lensi ng density ( Eg. |3l) 



In Coe et al. 



pOOSI), we tested our method's ability 
to recover a mass distribution similar to A1689 given 
93 multiple images of 19 lensed galaxies. Our recovered 
mass profile matched the input mass profile extremely 
well. In future work (§7|, we will quantify our ability to 
recover input mass profiles and concentrations. 

In Fig. [9] (left), we verify that NFW fits to our M{< R) 
are more tightly constrained than NFW fits to our k{R). 
In the right panel we provide an illustrative explanation. 
There is a wider range of NFW profiles which fit k{R) 
well (x^ < 1) than fit M{< R) well0 

Note that mass profiles of simulated galaxy clusters 
are often quoted in terms of their 3-D mass density p{r). 
Lensing alone cannot measure p{r) (althoug h the ad 



dition of other data can constrain p(r), as in Morandi 
et al.|2010J . Reported measurements of M(< A) m sim- 
ulations would enable more direct comparison between 
lensing observations and simulations 



Along those lines, Broadhurst & Barkana (2008) com- 
pared observed measurements of Einstein radii Re to 
those measured in simulations. The Einstein radius is 
a robust quantity which, for an axisymmctric potential, 
can be measured roughly by eye from the positions (and 
redshifts) of large arcs, then refined by model fitting. 
Measurement of Re yields the equally robust M(< Re), 
as the average surface density within Re must be equal 
to the critical lensing density T,crit (as a function of red- 
shift). Our shift from analyzing k(R) to the more ro- 
bust M{< R) measurements (adopted by other authors 
as well) is a natu ral extension of the shift to ro bust Re 
measurements by Broadhurst & Barkana ( 2008 ) . 



6.2. Simultaneous fits to the SL and WL data 

Lensing-based constraints on mass concentration are 
best derived from simultaneous SL + WL analyses. To- 
gether, SL and WL probe a sufficient range of radius from 
the cluster center, capturing the profile turnover from 
p (X r~^ and M(< R) oc R^ near the core to p ex r~^ an d 



M(< R) ^ rising slowly in the outskirts (see Coe||2010" I 
Meneghetti et al. (20091) quantify these statements, show- 



ing that KL-I-WL analyses of simulated clusters yield con- 
centration measurements to ~ 3x greater precision than 
WL-only analyses (11% scatter versus 33% scatter) and 
^ Sxgreater precision than SL-only analyses (59% scat- 
ter) ^ 

In i'ig. [To] we show NFW and Sersic profiles fit si- 
multaneously to SL M(< R) from our model and WL 
shears measured in previous works (see below). We ob- 

j+54 



258^^* kpc, C200 = 9.2 ± 1.2 



tain an NFW fit of r, 
(cvir — ^^-^-ia) and a Sersic fit of Re — 273.6 kpc, 
Ke = 0.2939, n = 1.425 (6„ = 2.524). Also plotted are 
fits to the SL data only. 



= 1.8 



lO^^Mcr^hZ 



Our NFW fits yield M200 ~ -^-"-os ^ ^^ -'"©"-to 
within r2oo = ^A^qI Mpc h^j^ and Af„,> = 2.0+o;3 x 
IO^^MqIi^q within the virial radius ryir = 3.0 ± 
0.2 Mpc hjQ. By definition, the average overdensities 
within these regions are Ac = 200 and 115, respectively, 
relative to the critical density to close the universe (see 
Appendix [a]). 

We fit to an ensemble of WL data measured in ACS 
hereafter M07 



20071 



Leonard et al. 



(Mede zinski et al 

.2007, hereafter Le(J7), Subaru (.Broadhurst et al.||2005b 



Hereafter B05b; M07), and CI-'HT images (L07). We 
also experiment by fitting to various subsets of this data 
ensemble. 
Fig 
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compares constraints on {Re, C200) for NEW 
fits to'HL and WL combined using three different WL 
data sets. We consider the L07 CFHT data (right), this 
data combined with the M07 ACS data (middle) , and all 



12 is similar but shows constraints on 



data (left). Fig. 

{M2QO, C20o)- 

Our derived SL+WL parameters are not affected 
greatly by our choice of WL data subset. However we 
note that inclusion of the Subaru WL measurements does 
increase the concentration slightly. 

Note that B05b converted their shear measurements 
7(i?) to mass measurements k{R) using their measure- 
ments of magnification ij.{R) to break the mass-sheet de- 
generacy. For the purposes of fitting an NFW (or Sersic) 
profile to the data, we prefer to avoid this step which 
may introduce additional uncertainty. HOG use the more 
direct approach which we use here as well. 

6.3. Comparison to Previous Work 

In Fig.ITSJwe compare our NFW fit parameters to pub- 
lished values from other similar studies: fits to SL, WL, 



'^^ One might worry about correlated uncertainties in our Af (< 
R) bins. This would be a concern only if k{R) were a more funda- 
mental observable than M(< R). We believe the opposite to be the 
case; lensing most fundamentally constrains M(< R). One should 
instead worry that derived measurements of k{9) are correlated; 
they must add to produce the correct M(< R) to deflect images to 
their observed radii. 



^^ Meneg hetti et ah] (| 2009) also find that concentrations derived 
from yL-only analyses are biased high by ~ 60% on average. They 
attribute this bias to their multi-component fit. They include 
an isothermal BCG component, but the BCG profile is actually 
steeper, and the parent halo must compensate with a higher con- 
centration. Care must be taken in these comparisons however, as 
we and others measure concentration by fitting to the total mass 
profile (parent halo plus galaxy halos), not the DM halo alone. 
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Fig. 9. — Fits to the enclosed mass Af(< R) provide tighter constraints than fits to the mass surface density k{R). Left: NFW fit 
confidence contours (l-cr and 2-cr) in the plane of Einstein radius Re{zs = 2) and concentration C2oo- Best fit values are plotted as X's. 
(Note the NFW scale radius Ts is fixed by each (-Rg, c) pair.) Right: Illustration of why M{< R) provides tighter constraints than k{R). 
Right Top: Best NFW fit (green) to the observed k(R) and two other fits with reduced x^ ~ 1. While these are all reasonable fits to k(R), 
the red (blue) curve significantly overestimates (underestimates) the observed M{< R) (grey; right bottom). 
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Fig. 10. — NFW and Sersic profiles fit simultaneously to SL and WL data. Also plotted are fits to the SL data alone. Left: The SL mass 
profile {R < 70") derived in this work. Right: WL reduced tangential shears (R > 70") measured in ACS (M07, Lc07), Subaru (B05b, 
M07), and CFHT images (L07). 



or both simultaneously. These published values and more 
can also be found in Tabled 

Our concentration derived from SL alone (c2oo — 7.6 ± 
1.3) is higher than that obtained in previous SL analyses 
(c2oo ~ 6), though our values agree roughly within the 1- 
a uncertainties. In Fig. [MJ we comp are ou r SL A/(< R) 
to that obtained by [HalEola et aL (2006). Just inside 
the Einstein radius, our M{< R) is steeper and claims a 
lower uncertainty. This budges the concentration value 
higher. 

We have neglected to mask out substructure which can 
potentially be problematic for NFW fits (e.g., Hennawi 
et al.|[2007 ). However we note the main subcluinp which 
we might exclude is located between ^ 150 - 200 kpc 
from the cluster center. Our mass excess relative to H06 
is at smaller radius, between ^ 100 - 150 kpc. 

From our SL+WL fits, we find Caoo = 9.2 ± 1.2. This 
lies between the values of C200 — 7.6^0 5 found by H06 

and C200 = lO-S^o's found by B05b. We note our re- 
sults are in good agreement with more recent analyses 



that have found v alues of 0200 R^ 10 ± 0.7 from SL+WL 



number counts (Umetsu & Broadhurst 2008[ Umetsu 
et al . 1 120091) and SL+WL + X-ray data (ILenize et al. 

m^ — 

We note our uncertainty contours closely follow the 
Re = 47"5 isocontour of constant Einstein radius. The 
best fit NFW parameters obtained in previous SL or 
SL+WL studies all fall along this line, as the Einstein 
radius is tightly constrained by SL. 

6.4. Tension between the SL and WL data? 

In every SL+WL study pubhshed to date for A1689, 
including ours, the WL data prefer a higher conccntra- 
tion and steeper mass slope than the SL data (Figs. 



and 15). While we have found a single NFW profile 
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which provides a decent fit to both the SL and WL data, 
there is a slight tension between the two. While com- 
patible at l-cr depending on the WL data set chosen (see 
Fig. 11), the WL signal is a bit too low and/or falls off 

Alternatively, 
We might 



a bit too quickly relative to the SL data, 
the inner SL mass profile is a bit too shallow. 
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Fig. 11. — Constraints on Einstein radius Re and concentration C200 from an NFW profile fit to the observed SL (red), WL (green), and 
both combined (yellow). Three sets of WL data are used - left: all; middle: CFHT (L07) & ACS (M07); right: CFHT (L07). Best fits are 
marked with X's. Confidence contours are l-c (68%, Ax^ = 2.3) and 2-a (95%, Ax^ = 6.17). A black dashed line also indicates the 1-D 
1-(T WL contours (Ax'^ = 1), which give the 68% bounds if one parameter is marginalized over. Isocontours of M200 S'l'e overplotted. 
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Fig. 12. — Same as Fig. |ll| but for mass within the virial radius M200 and NFW concentration C200- Isocontours of Re (zs 
are overplotted. 



: 2; arcsec) 



suspect the WL signal has been diluted from contami- 
nation of the sample by unlensed (foreground or cluster) 
objects. However M07 carefully considered and quanti- 
fied dilution in their work. 

We note that these deviations from NF W are qualita- 
tively in ac cord with those proposed by |Lapi fc Cava- 
liere (2009), as found in their semi-analytic modeling of 
cluster formation. This could prove interesting if similar 
variations are detected in many clusters. Slight devia- 
tions from NFW are not surprising in a single cluster, as 
simulated cluster profiles do exhibit intrinsic variations. 

We note that L07 had claimed to resolve the discrep- 
ancy of SL versus WL concentrations. They fit NFW 
profiles separately to their SL and WL data, finding 
concentration parameters of C200 = 6.0 ± 0.6 (3-cr) and 
C200 = 7.6 ± 1.6 (l-cr), respectively. As the error bars 
overlap, they claimed agreement. However as we showed 
in Fig. [11] SL-only and WL-only NFW fits may yield 
similar concentration values but very different Einstein 
radii. L07's best fit to the WL data yields Re ~ 30" 
{zg = 2), significantly lower than the value Re ~ 47" 
tightly constrained by SL. This fit unde rpre dicts mass at 
all radii within the Einstein radius (Fig. 13 ). L07 did not 
attempt to provide a single profile which provides an ac- 
ceptable fit to both the SL and WL data simultaneously. 

6.5. Comparison to Predictions 

Our results support findings that A1689 has a higher 
concentration than predict ed by simulations fo r a cluster 
of its mass. According to Duffy et al. (2008), a cluster 
such as A1689 with M200 w 2 x lO^'^M© 



should have C200 = 3.0_Q;g {cyir = 3.9_ii). Relaxed 
clusters are more symmetric yielding better fits to NFW 
profiles and ~ 20% higher concentrations for clusters of 



this mass: C200 
from 



3.4 



Maccio et al. (2008) 



+1.4 
-1.0 



(c. 



= 4.4 



+1.8 



vir — ^-^-1.3, 

similar (~ 



are 

v^-l 



. Estimates 
10% lower) 



though the normal c oTJT+z) ^ scaling relation (Bullock 
et al ||2001[) drags the predictions down another ~ 10%. 
Duffyetal.| ( |2008[ ) instead find roughly C200 oc {l+z)-°-^^ 

and Cvir oc (1 -|- z)~ °'^°. 

Pred ictio ns from IBuUock et al. ( 2001 ) 
( |2007t , and |Netoet al ( 2007^ 



Gentile et al. 
wTio analyzed the Millen- 



Duffy et al. 
1. 'i'he Mil- 



nium simulation) are also fairly similar to the 
(|2008| predictions despite concerns about eac 
lenmum simulation used a WM APl cosmolo g y ([Spergel 
et al.p2003|) in cluding as 



0.9. Duffy et al. (20DSraHd 



■1 

m 



Maccio et al. (2008) found this res ults m concentrations 
15% higher tlTairtheir WMAPS QKomatsu et al.|2009 l 
input erg = 0.796 p] The|Bullock et al.|(|2001[) simulations 



used gs = 1-0 ana did not produce halos as massive as 



A1689. Zhao et al. ( 2003 ) cites the dangers of extrapolat- 



ing these results to higher mass. Gentile et al. (2007| is 
a reformulation of the original Navarro et al. (1996) pre- 
scription to WMAP3 ( ^Spergel et al.,2007) . Those early 
simulations were lower resolution and prod uced only 19 
halos (co mpared to ^1,000 a nd ^10,000 for|Duffy et al. 
2 008 1 and [Maccio et al l 20081 respectively) 

— mit — 



Hennawi et al. 



at z 



0.186 



measure significantly larger con- 
centrations in their simulations. For A1689, they predict 

^® This value is in excellent agreement with the WMAP 7-vear 
maximum likelihood value erg = 0.803 (Komatsu et al.|2010||. 
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Fig. 13. — NFW fit parameters (r^, C200 [cuir]) found in tfiis 
paper and in previous works compared to those predicted from 
simulations. Filled circles give best fit SL parameters from BOS, 
Z06, H06, L07, and ours Vs = 338^7°^ kpc, C200 = 7.6 ±1.3 (c„i,. = 
9.6+1-^ 
and ours r. 



■ 33St\T kpc, C200 
-^ g). Filled upright triangles give SL+WL fits from B05b, H06, 



:258"' 



For 



»t4o kpc, C200 = 9-2 ± 1.2 {c^,r = ll-5ll;i) 
clarity, error bars are plotted vertically, though they should actu- 
ally all tightly follow the Re = 47'.' 5 {zg = 2) isocontour, as shown 
for our SL+WL fit. We also plot the WL fit obtained by L07. This 
fit, with Re = 30" (zs = 2), is the only fit to fall far from the 
Re = 47'.' 5 (23 = 2) isocontour. (Again note the error bar is not 
properly oriented.) Along the isocontour M200 = 2 X lO^^Af© (our 
best fit value), we have plotted concentrations observ ed in simula- 
tionsBuUock et al. 2001 Hcnnawi et al. 2007 Neto et al-||2007 



2001 Hcnnawi et al. 2007 ii__i 

|Maccio e t al. 20U8, Uutty et al. 20U8), mcludmg the origmal INl' W 
pr escrip tion (Navarro et al. 1996) adapted to the present cosmology 



r escrip l 
AGent 



bylGentilcct al. (2007). All assume a scaling of c oc (1 + 2:) ex- 
cept [Duffy "et al. (2008) which find roughly c oc (I + 2) °'*5. Their 
shallower dependence on z happens to bring their lower values back 
in line with some earlier predictions. Filled symbols are plotted for 
relaxed subsets of the full samples. Again the error bars are plotted 
vertically, but here should follow the M200 = 2 X lO^^Af isocon- 
tour. Finally, the mass-concent ration relation observed in clusters 
||Comerford &: Natarajan|2007[) is plotted as the yellow square. 
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Fig. 14. — Comparison of our SL model M(< R) to that of 
H06. Also plotted are four different NFW profiles with C200 = 
6,7,8,9, all with the correct Einstein radius. All of the profiles are 
normalized to the C200 = 6 NFW profile. Our SL model has a 
higher concentration (c200 = 7.6 it 1.3) than that of H06 (c200 = 
6 ±0.5). 



C200 = 4.4 



Q^, ^ 50% larger than predicted by Duffy 
(2008|). Their use of erg = 0.95 probably oiily re 

~ 20%. The remaining 



et al. 

suits m concentrations inflated by 

disagreement may be a result of their halo density fitting 

procedure which they claim is better for comparison with 

lensing measurements. 




100 150 

Radius (l<pc) 

Fig. 15. — Our SL-derived mass density profile k{R) (black 
squares), along with our NFW profile fits and those published else- 
where. All fits which incorporate SL data (dashed: SL only; solid: 
SL+WL) provide reasonable fits to our SL data at most radii. The 
L07 WL-only NFW fit (dotted line) underpredicts the mass surface 
density at all radii plotted here. 



All of these predictions (detailed further in Coe||2010" I 



are plotted in Fig. [13 The predictions are all signih 
cantly lower than oursL-|-WL-derived C200 = 9.2 ± 1.2. 
Clusters selected with a lensing bias may h ave mea- 
sured conce ntratio ns bi ased high by ^ 34% Hennaw 
eFaI| ( [2007|), ^ 50% [OK uri fc Blandford. ( .2009j, or even - 
TUD% [JVlenegh etti et al.| ( I201U| ), on average. I'his results 
from a combination ot higher intrinsic (3-D) concentra- 
tions and additional mass along the line of sight (either 
due to cluster elongation or otherwise) resulting in higher 
projected (2-D) concentrations. The C200 = 4.4lj!)g pre- 

100% bias might 
9. However such 



dieted by Hennawi et al. ( 2007 ) plus a 
begin to explain measurements of C200 
extremely high biases are expected only for less massive 
clusters. Even accounting for such biases, the high con- 
centration of A1689 seems unlikely given the results from 
ACDM simulation s ( [Broadhurst fc Barkana||2008[ [Oguri 
fc Blandford|[2009| . 

i'here are fewer ]Dublished fits of S ersic profiles to simu- 



lated cluster halos. We do note that Merritt et al. (2005 1 



found n — 2.38±0.25 for their cluster sample]^ ^fOur best 
fit 71 = 1.425 is much lower giving our density profile a 
more rapidly varying slope. This allows it to fit both 
the SL profile and the low WL signal measured by some 
authors at large radius. We have not experimented with 
fitting Sersic profiles to various subsets of the WL data. 
Recently the Einast o profile has gained popularity over 
the Sersic profile (e.g., Navarro et al. 2010 ). The two have 
similar forms, but the former gives density p(r) as a func- 
tion of 3-D radius, while the latter gives surface density 
I](i?) as a function of projected 2-D radius. The former 
was found (Merritt et al. 2006) to provide better fits to 
halos of a wide range of masses, though the latter per- 
formed slightly better specifically for cluster halos (and 
only slightly worse for galaxy halos). Mass concentration 
relations derived from Einas t o profile fits have b een pub- 
lished by [Duffy et al.| |2008[ ); |Gao et al.| ( |2008[ ) ; [Hayashi 



3^^ Alternativel y, fitting a dcproiectcd Sersic profile to the 3D 
^tial density p( r),|Merritr et al. (2005) and Merritt et al. (200^ 
Found riri = 2.99 ± 0.49 and n^ = 2.89 ± 0.49, respectively. Ttie 
deprojcction tfiey used (Prugnicl & Siniicn 1997) is approximate 
land thus does not yield exactly the same n as fitting directly to 



[the surface density. 
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fc White] p)08| ). These fits yield sli ghtly (< 2 0%) dif- 
terent concentrations than NFW fits ( Coe||2010 |. We do 
not explore Einasto profile fits here. 

The Sersic profile remains int riguing for the di rect com- 
parisons which can be made ( jMerritt et al. 2005 ) to many 
pubhshed Sersic fits to galaxy lumin osity profiles, (al- 
thoug h this may just be coincidental IDhar fc Williams| 
20101. We also note the possibilit y of using well est ab- 
UsHed software such as GALFIT ( [Peng et aL]|2002[ ) to 
derive Sersic parameters for surface density mass maps. 

7. FUTURE WORK 

In this paper, we have focused on measurement of the 
radial mass profile of A1689. However, a key strength 
of LensPerfect is its ability to map massive substructure 
without assuming LTM. In future work, we will verify 
in detail our ability to resolve halo subclumps as well as 
measure their masses. The latter proves difficult (both 
in observed and simulated halos) as the subhalo masses 
must be disent angled cleanly from the g reater parent halo 
(e.g., [Natarajan et al.||2Q09( [Julio fc Kneib 2009J. 



One novel method developed recently by [Powell et al.| 
( 2009 1 demonstrates the ability to detect subclumps in a 
2-D mass map down to 10^"^ or even 10^^ Mq and mea- 
sure their masses to within a factor of two. By weighing 
subhalos associated with cluster galaxies, we may provide 
evidence for galaxy halo stripping in individual galaxies. 
This would provide an excellent complement to studies 
which have measure d stripping "globally" , averaged over 
many galaxy halos (Natarajan et al.'1998' '2002'; Gavazzi 
et_ all 12004; Limousin et al. 2007; Halkola et al. 200"" 



er 



Natarajan et al. 2009). 

~ |WiUiams fc Saha, ( ,2004| and |Sahaet aL|(|2007[ ) have ex- 
p erimented with several met hods to identify subst ruct ure 
in their non-LTM (PixeLens) mass models. From the 2-D 

mass density map k(6'), they have subtracted each of the 
following: the average k{R) in that radial bin; n{6) 180° 
across (directly oppos ite w.r.t. the cente r); or the best fit- 
ting NFW profile. In|Saha et al.| (|20q7|), they show that 
the observed substructure (extcndedrneso-structure" ) 
appears to correlate with the luminous galaxies. 

In future work we will develop algorithms to thor- 
oughly explore the range of model solutions which per- 
fectly reproduce all observed multiple image positions. 
This method must take care not to remain trapped in a 
local minimum near our "best" solution. Also, we must 
correctly account for the larger uncertainties in voids be- 
twee n the multiple im ages. 

(2008), we discussed our ability to mold 

by adding extra artificial constraints. 

mass models would successfully repro- 



al. 



In Coe et 
the mass map 
These modified 
duce all of the input data plus the artificial constraints. 
These added constraints squeeze the mass model, tweak- 
ing the positions of subclumps or increasing the concen- 
tration of mass in the desired regions. We must explore 
such solutions to accurately account for all uncertainties. 
We may also mold our mass models in attempts to force 
mass to follow light more closely. We might construct 
that mass model which follow s light best, as in (for ex- 
ample jSaha fc Williams||1997[ ) . 

Ultimately a hybrid approach combining non-LTM and 
flexible-LTM components may prove ideal. A prior of 
LATM (light approximately traces mass) could be as- 



sumed. This ideal method would include a parent halo, 
galaxy components, and line of sight structure, all with 
sufficient flexibility. Each galaxy component might be 
allowed to vary individually in M/L, radial scale, trun- 
cation radius, and perhaps position. Different forms may 
even be explored: truncated isothermal ellipsoid versus 
NFW, for example. The parent halo should be very flexi- 
ble (a multi-scale grid perhaps) to allow for the asymme- 
tries induced as galaxies infall and their stripped mass is 
strewn about the cluster. Additional mass planes behind 
the lens should also be modeled. 

Future deeper observations of galaxy clusters such as 
A1689 may reveal hundreds of multiple images. This 
wealth of constraints will allow truly high-definition mass 
models, which clearly resolve individual galaxy halos and 
perhaps dark subhalos as well ( ICoc 20 09 ) . As more mul- 
tiple images are revealed and greater details are obtained, 
line of sight structure will need to be taken into account. 
While seemingly a nuisance, this raises the prospect of 
strong lensing tomography - the mapping of mass in mul- 
tiple lens planes. 

Over the next 2 to 3 years we are looking forward to 
new HST images from the CLASH Multi-Cycle Treasury 
Program. As mentioned in the introduction, this pro- 
gram will image 25 clusters to a depth of 20 orbits each, 
equal to the depth of the A1689 images studied here. 
Analysis of these ACS and WFC3 images along with 
supporting data will yield measurements of mass con- 
centration for a sizable sample of clusters, selected free 
of lensing bias. By comparing these values to the con- 
centrations of simulated halos, we expect to either show 
agreement or detect average deviations as small as 15% 
with 99% confidence. 

8. SUMMARY 

We have presented a strong lensing mass model of 
A1689 which resolves structures down to ^ 25 kpc on 
average within the central 400 kpc diameter without as- 
suming light traces mass. The most luminous galaxies 
appear to trace the mass distribution fairly well but with 
some deviations which may prove interesting, pending 
verification. This is the highest resolution mass map of 
any galaxy cluster to date. Our mass model perfectly re- 
produces the observed positions of 168 multiple images 
of 55 knots within 135 images of 42 galaxies. Included 
are 20 new candidate multiple images of 8 galaxies which 
we have identified in this work. We have also tweaked 
some of the identifications from previous works, discard- 
ing three suspect central images and adding one new one. 

Compiling published weak lensing measurements from 
ACS, Subaru, and CFHT images, we find that a single 
mass profile, either NFW or Sersic, is able to provide a 
decent fit simultaneously to both the observed weak and 
strong lensing. However there remains a slight tension in 
that the weak lensing data prefer higher concentrations 
than the strong lensing data. 

Based on simultaneous fitting of the strong and weak 
lensing data, we measure a NFW central mass concen- 
tration of C200 = 9-2 ± 1.2. Thus we concur with previ- 
ous claims that the mass profile of A1689 appears to be 
more centrally concentrated than clusters of similar mass 
(Af2oo -- 2 X IQi^Mo/ifo^ = 1.4 x W^M^hr'^) formed in 
CDM simulations (c2oo ^ 3). 

In future work, we will perform further simulations and 
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verify our ability to not only resolve substructure but 
measure the masses of these halo subclumps. Backed 
by these tests, future analyses of our substructure maps 
should yield further evidence for galaxy halo stripping 
in cluster environments, test the degree to which light 
traces mass, and perhaps even lead to detection of dark 
subhalos, should they exist. 
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TABLE 4 

Spectroscopic Redshifts of galaxies within the ACS field of view of 

Abell 1689 used for our photometric redshift tests 



ID 


RA & DEC (J2000)=' 


Survey'' 


^spec 


451 


13:11:32.606 


-1:19:22.920 


Miekse 


0.0000 


366 


13:11:24.802 


-1:20:23.070 


Frye 


0.0000 


374 


13:11:26.919 


-1:20:00.650 


Frye 


0.0000 


216 


13:11:28.344 


-1:18:32.380 


MOS 


0.0130 


213 


13:11:27.224 


-1:20:10.180 


MOS 


0.0862 


220 


13:11:33.034 


-1:21:24.780 


MOS 


0.1430 


74 


13:11:25.382 


-1:20:17.140 


Balogh* 


0.1660 


7 


13:11:32.089 


-1:19:36.600 


Balogh* 


0.1690 


127 


13:11:28.194 


-1:20:43.680 


MOS* 


0.1712 


194 


13:11:30.934 


-1:20:30.680 


Teague* 


0.1739 


387 


13:11:30.508 


-1:19:34.670 


Frye 


0.1740 


378 


13:11:27.851 


-1:20:07.650 


Frye 


0.1750 


192 


13:11:30.234 


-1:20:27.680 


Teague* 


0.1750 


124 


13:11:27.894 


-1:21:36.880 


MOS* 


0.1752 


184 


13:11:26.884 


-1:19:37.180 


Teague* 


0.1754 


43 


13:11:37.981 


-1:20:09.780 


Balogh* 


0.1760 


108 


13:11:24.284 


-1:21:14.480 


MOS* 


0.1766 


198 


13:11:32.094 


-1:21:38.380 


Teague* 


0.1770 


4 


13:11:28.344 


-1:19:58.300 


Balogh* 


0.1775 


154 


13:11:32.464 


-1:22:18.280 


MOS* 


0.1794 


112 


13:11:25.104 


-1:19:31.380 


MOS* 


0.1797 


383 


13:11:28.886 


-1:20:01.940 


Frye 


0.1800 


141 


13:11:29.994 


-1:20:17.580 


MOS* 


0.1801 


152 


13:11:32.104 


-1:19:47.180 


MOS* 


0.1801 


61 


13:11:31.510 


-1:19:24.870 


Balogh* 


0.1817 


9 


13:11:34.058 


-1:21:01.990 


Balogh* 


0.1826 


370 


13:11:26.237 


-1:19:56.450 


Frye 


0.1830 


201 


13:11:35.374 


-1:20:43.080 


Teague* 


0.1835 


64 


13:11:29.386 


-1:18:34.790 


Balogh* 


0.1836 


172 


13:11:37.834 


-1:19:20.880 


MOS* 


0.1839 


375 


13:11:27.148 


-1:18:48.440 


Frye 


0.1840 


135 


13:11:29.474 


-1:20:28.080 


MOS* 


0.1842 


118 


13:11:27.084 


-1:18:48.880 


MOS* 


0.1852 


153 


13:11:32.204 


-1:22:10.980 


MOS* 


0.1855 


188 


13:11:29.074 


-1:21:37.380 


Teague* 


0.1858 


132 


13:11:28.724 


-1:19:02.980 


MOS* 


0.1859 


450 


13:11:32.761 


-1:19:48.930 


Miekse 


0.1859 


63 


13:11:29.869 


-1:20:15.230 


Balogh* 


0.1868 


165 


13:11:35.344 


-1:21:33.480 


MOS* 


0.1870 


69 


13:11:28.023 


-1:18:43.890 


Balogh* 


0.1870 


110 


13:11:24.434 


-1:21:11.180 


MOS* 


0.1870 


147 


13:11:31.064 


-1:21:27.980 


MOS* 


0.1872 


384 


13:11:29.100 


-1:19:46.920 


Frye 


0.1880 


6 


13:11:31.400 


-1:19:32.840 


Balogh* 


0.1884 


196 


13:11:31.074 


-1:20:52.780 


Teague* 


0.1885 


419 


13:11:38.051 


-1:19:58.250 


Frye 


0.1890 


390 


13:11:31.472 


-1:21:05.940 


Frye 


0.1890 


163 


13:11:34.754 


-1:20:59.480 


MOS* 


0.1895 


189 


13:11:29.114 


-1:21:55.480 


Teague* 


0.1908 


131 


13:11:28.594 


-1:20:26.780 


MOS* 


0.1909 


70 


13:11:26.812 


-1:19:43.090 


Balogh* 


0.1910 


146 


13:11:30.564 


-1:20:43.980 


MOS* 


0.1918 


1 


13:11:25.354 


-1:20:37.060 


Balogh* 


0.1922 


133 


13:11:29.294 


-1:19:16.980 


MOS* 


0.1932 


187 


13:11:28.994 


-1:21:16.980 


Teague* 


0.1947 


71 


13:11:25.925 


-1:19:51.950 


Balogh* 


0.1950 


119 


13:11:27.094 


-1:21:43.080 


MOS* 


0.1955 


12 


13:11:36.616 


-1:19:42.800 


Balogh* 


0.1960 


148 


13:11:31.124 


-1:21:25.380 


MOS* 


0.1972 


123 


13:11:27.834 


-1:21:13.080 


MOS* 


0.1977 


115 


13:11:25.974 


-1:19:35.680 


MOS* 


0.1983 
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RA & DEC (J2000)'' 



Survey ° 



140 


13:11:29.984 


-1:22:07.580 


MOS* 


0.1985 


11 


13:11:35.595 


-1:20:12.400 


Balogh* 


0.1995 


128 


13:11:28.324 


-1:18:45.180 


MOS* 


0.1999 


395 


13:11:33.234 


-1:19:17.000 


Fryc 


0.2000 


139 


13:11:29.974 


-1:20:40.480 


MOS* 


0.2002 


158 


13:11:32.704 


-1:19:32.080 


MOS* 


0.2009 


143 


13:11:30.044 


-1:20:42.980 


MOS* 


0.2012 


452 


13:11:29.905 


-1:20:05.430 


Miekse 


0.2014 


157 


13:11:32.654 


-1:19:58.980 


MOS* 


0.2022 


2 


13:11:27.109 


-1:20:58.420 


Balogh* 


0.2147 


82 


13:11:23.065 


-1:21:17.280 


Balogh* 


0.2150 


3 


13:11:27.681 


-1:21:07.190 


Balogh* 


0.2158 


398 


13:11:33.555 


-1:19:01.500 


Frye 


0.2440 


221 


13:11:33.724 


-1:19:39.180 


MOS 


0.3100 


214 


13:11:27.844 


-1:18:52.980 


MOS 


0.3840 


399 


13:11:33.621 


-1:22:02.190 


Frye 


0.3870 


376 


13:11:27.204 


-1:18:49.920 


Frye 


0.4800 


365 


13:11:24.652 


-1:20:03.380 


Frye 


0.4810 


405 


13:11:35.088 


-1:21:26.010 


Frye 


0.5840 


406 


13:11:35.258 


-1:20:30.170 


Frye 


0.5870 


418 


13:11:37.699 


-1:19:50.030 


Frye 


0.6250 


388 


13:11:30.677 


-1:18:55.500 


Frye 


0.6760 


402 


13:11:34.328 


-1:19:05.100 


Frye 


0.6760 


389 


13:11:30.751 


-1:21:38.780 


Frye 


0.6910 


379 


13:11:28.221 


-1:20:50.930 


Frye 


0.7030 


380 


13:11:28.325 


-1:18:27.500 


Frye 


0.7100 


367 


13:11:24.960 


-1:19:36.610 


Frye 


0.7220 


364 


13:11:24.609 


-1:19:20.830 


Frye 


0.7570 


411 


13:11:36.536 


-1:19:25.030 


Frye 


0.7900 


396 


13:11:33.028 


-1:19:14.650 


Frye 


0.7900 


217 


13:11:28.764 


-1:21:43.580 


MOS 


0.7900 


414 


13:11:37.183 


-1:20:16.890 


Frye 


0.8130 


362 


13:11:24.240 


-1:19:52.680 


Frye 


0.8570 


363 


13:11:24.398 


-1:19:36.880 


Frye 


0.8950 


409 


13:11:36.052 


-1:19:24.710 


Frye 


0.9160 


413 


13:11:37.089 


-1:19:26.100 


Frye 


0.9240 


373 


13:11:26.683 


-1:19:36.770 


Frye 


0.9590 


393 


13:11:32.056 


-1:21:55.410 


Frye 


0.9590 


356 


13:11:22.472 


-1:20:39.220 


Frye 


0.9600 


401 


13:11:34.242 


-1:19:23.970 


Frye 


1.0510 


377 


13:11:27.191 


-1:18:26.550 


Frye 


1.1120 


361 


13:11:24.186 


-1:19:56.550 


Frye 


1.1550 


400 


13:11:33.950 


-1:19:15.750 


Frye 


1.3620 


306 


13:11:33.981 


-1:20:50.950 


Broadhurst 


1.3700 


302 


13:11:27.359 


-1:20:54.900 


Broadhurst 


1.8200 


301 


13:11:30.294 


-1:19:50.920 


Broadhurst 


1.8200 


371 


13:11:26.541 


-1:19:55.700 


Frye 


2.5360 


307 


13:11:33.065 


-1:20:27.450 


Broadhurst 


3.0410 


308 


13:11:26.441 


-1:19:56.940 


Broadhurst 


3.0410 


385 


13:11:29.969 


-1:19:14.880 


Frye 


3.8500 


368 


13:11:25.447 


-1:20:51.740 


Frye 


4.8680 


404 


13:11:34.994 


-1:19:51.170 


Frye 


5.1200 



Slight modifications have been made to t he prev iously published coordinates where ne c essarv- 
Redshift surveys are 



Balogh et al. 



I 2002 I, MOS I Due et al.|2002 i 



Teague et al, 



11990 



Broadhurst et al. 



I 2005a I, and 



Frye et al 


2002 


2007| 



Asterisks (*) indicate galaxies designated as "spectroscopically-confirmed cluster members" in each survey. 
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TABLE 5 
Multiple Images Produced by Abell 1689 



ID 



10=* 
BOS 



X" 

(pixel) 



(pixel) 



R.A. 
(J2000.0) 



Decl. 
(J2000.0) 



^ 775 
(AB mag) 



Ref.' 



Zb- 



2 f 
■^mod 



10 



11 



12 



13 



14 



15 



16 



17 



a 


2213 


1055 


13 11 26.452 


-1 19 56.75 


23.309 ± 0.008 


3.04 


2,3 3.2i+;;:5,' 


b 


2130 


1040 


13 11 26.289 


-1 20 00.19 


23.595 ± 0.009 


(3.04) 




3.88 ±0.48 


c 


1356 


2557 


13 11 29.773 


-1 21 07.43 


24.509 ±0.017 


3.04 


2,3 3.79 ±0.47 


d 


2499 


3112 


13 11 33.066 


-1 20 27.47 


24.033 ±0.011 


(3.04) 




3.76 ±0.47 


e 


2745 


2621 


13 11 31.932 


-1 20 05.91 


24.638 ±0.018 


(3.04) 




3.i7t«:I? 


f 


1890 


2333 


13 11 29.852 


-1 20 38.50 


24.942 ± 0.050 


(3.04) 




o.5it^:«? 


a 


2245 


1064 


13 11 26.524 


-1 19 55.49 


23.238 ±0.013 


2.53 


3 


16+"™ 


b 


2522 


3069 


13 11 32.969 


-1 20 25.51 


23.912 ±0.015 


(2.53) 




0.20l°f2 


c 


2728 


2644 


13 11 31.978 


-1 20 07.17 


24.285 ± 0.024 


(2.53) 




2m+_ffs 


d 


1386 


2556 


13 11 29.812 


-1 21 06.05 


24.351 ± 0.025 


2.5 


4 


2.661^4^ 


e 


1876 


2349 


13 11 29.881 


-1 20 39.48 


24.806 ± 0.041 


(2.53) 




03ltli\ 


a 


2372 


2832 


13 11 32.041 


-1 20 27.27 


26.099 ± 0.056 






0.96 ±0.19 


b 


2279 


2921 


13 11 32.178 


-1 20 33.37 


26.287 ±0.054 






5.22 ±0.61 


c 


1809 


2984 


13 11 31.703 


-1 20 55.99 


27.243 ±0.129 






0.81+^-p 


a 


1844 


3124 


13 11 32.175 


-1 20 57.37 


24.172 ±0.017 


1.1 


4 


o.79t°;t^ 


b 


1369 


2801 


13 11 30.528 


-1 21 12.02 


23.534 ±0.015 


(1.1) 




1 r|n + 0.53 


c 


2553 


2322 


13 11 30.758 


-1 20 08.25 


25.021 ±0.027 


(1.1) 




1.061^:^0 


d 


1492 


1338 


13 11 26.285 


-1 20 35.40 


24.610 ±0.020 


(1.1) 




0.78i°:5i 


e 


2052 


2252 


13 11 29.837 


-1 20 29.45 


27.171 ±0.325 


(1.1) 




o-68i^i^ 


a 


1606 


2205 


13 11 29.064 


-1 20 48.64 


24.290 ±0.030 


2.6 


4 


2.911^4^ 


b 


1706 


2211 


13 11 29.224 


-1 20 44.24 


24.676 ±0.041 


(2.6) 




0-25tf,11 


c 


2751 


3343 


13 11 34.120 


-1 20 20.96 


25.103 ±0.027 


2.6 


4 


3.00 ±0.39 


a 


3097 


2066 


13 11 30.755 


-1 19 38.19 


23.521 ±0.015 


1.1 


4 


1.17±0.21 


b 


2811 


3058 


13 11 33.345 


-1 20 12.20 


23.624 ±0.016 


1.1 


4 


1.36 ±0.23 


c 


3055 


2744 


13 11 32.742 


-1 19 54.49 


22.131 ±0.006 


(1.1) 




0.98 ±0.19 


d 


2943 


2709 


13 11 32.478 


-1 19 58.81 


21.827 ±0.003 


(1.1) 




0.53±«t^ 


a 


1087 


1250 


13 11 25.446 


-1 20 51.87 


23.292 ±0.007 


4.87 


1 


4.93 ±0.58 


b 


2439 


2349 


13 11 30.678 


-1 20 13.99 


24.128 ±0.014 


4.8 


4 


4.86 ±0.57 


a 


1974 


3105 


13 11 32.302 


-1 20 51.09 


22.011 ±0.009 






n 1 5+0-31 


b 


1596 


2984 


13 11 31.402 


-1 21 05.63 


21.587 ±0.007 








U.48_o3g 


c 


2541 


2572 


13 11 31.495 


-1 20 14.10 


25.550 ±0.050 








0-22tg1^ 


d 


1674 


1001 


13 11 25.526 


-1 20 20.01 


22.278 ±0.005 








0.50 ±0.15 


a 


2850 


2032 


13 11 30.303 


-1 19 48.65 


25.661 ± 0.041 








4.80 ±0.57 


b 


2141 


3430 


13 11 33.519 


-1 20 50.42 


27.269 ±0.087 








0.79t^:?| 


c 


1072 


2347 


13 11 28.737 


-1 21 15.83 


25.560 ±0.036 








o.76tf,:r. 


d 


1645 


1264 


13 11 26.279 


-1 20 26.90 


26.941 ± 0.072 








4.95+2-g 


a 


2189 


3560 


13 11 33.980 


-1 20 51.01 


23.011 ±0.009 


1.83 


4 


0.26±^^? 


b 


2130 


1625 


13 11 28.055 


-1 20 12.61 


22.655 ±0.009 


(1.83) 




0.12±«:«' 


c 


2012 


2098 


13 11 29.316 


-1 20 27.99 


23.605 ±0.018 


(1.83) 




o.56±;-? 


a 


1824 


3522 


13 11 33.349 


-1 21 06.73 


23.794 ±0.015 


2.5 


4 


2.89 ±0.38 


b 


2462 


1801 


13 11 29.056 


-1 20 01.31 


23.402 ±0.015 


(2.5) 




2.821^:^^ 


c 


2062 


2135 


13 11 29.498 


-1 20 26.51 


23.729 ±0.015 


(2.5) 




0.53t°i° 


b 


1275 


1796 


13 11 27.361 


-1 20 54.94 


24.178 ±0.019 


1.83 


2,^ 


1 0.601^45 


c 


1311 


1730 


13 11 27.213 


-1 20 51.91 


23.845 ±0.015 


(1.83) 




0.60 ±0.16 


f 


1602 


1547 


13 11 27.072 


-1 20 34.86 


25.403 ± 0.075 


(1.83) 




9 in+o-48 

^■lU-0.60 


g • 


2197 


1563 


13 11 27.963 


-1 20 08.26 


22.833 ±0.005 


(1.83) 




-1 f,.+0.26 
l.b4_i 23 


h ■ 


2190 


3463 


13 11 33.688 


-1 20 48.91 


28.078 ±0.150 


(1.83) 




2.ootl:?« 


i 


2224 


1497 


13 11 27.802 


-1 20 05.64 


23.806 ±0.010 


(1.83) 




0-19±?,1^ 


a 


3610 


2512 


13 11 32.828 


-1 19 24.44 


23.632 ± 0.020 






1.03 ±0.20 


b 


3605 


2567 


13 11 32.986 


-1 19 25.83 


23.764 ±0.018 






0.37 ±0.13 


c 


3559 


2725 


13 11 33.398 


-1 19 31.27 


23.602 ±0.018 






1 ifi+0-33 


a 


639 


2648 


13 11 29.033 


-1 21 41.82 


25.141 ± 0.040 


3.4 


4 


n 55+2-66 


b 


677 


2772 


13 11 29.461 


-1 21 42.73 


25.700 ±0.052 


(3.4) 




o oy + 0.44 


a It 


.2 1356 


2557 


13 11 29.773 


-1 21 07.43 


25.508 ± 0.044 


1.8 


4 


1 91+0-45 


b IE 


.1 2213 


1055 


13 11 26.452 


-1 19 56.75 


24.722 ± 0.040 


(1.8) 




-, 70+0.84 


c 


2009 


2074 


13 11 29.239 


-1 20 27.62 


25.622 ±0.060 


(1.8) 




f. -0 + 1.96 


a 


1892 


1715 


13 11 27.990 


-1 20 25.29 


23.307 ±0.015 






l-89lo-28 


b 


1950 


1991 


13 11 28.905 


-1 20 28.53 


24.541 ± 0.049 






Qm+_lf, 


c 


2326 


3635 


13 11 34.400 


-1 20 46.40 


24.992 ± 0.031 






o.25lf,:r2 


a 


2240 


2437 


13 11 30.662 


-1 20 24.87 


24.133 ±0.021 


(2.6) 




r, 7r,+0.36 
^■'^-0.37 


b 


2152 


2389 


13 11 30.392 


-1 20 27.83 


23.077 ±0.010 


(2.e 


) 




40+0.14 



1.06 
1.08 
3.46 
2.37 
0.57 
2.22 
0.27 
0.47 
9.83 
0.90 
2.02 
4.30 
1.03 
3.72 
0.09 
1.91 
0.77 
0.80 
0.89 
0.03 
1.68 
1.18 
0.03 
0.73 
1.48 
4.20 
0.05 
1.75 
0.30 
1.36 
4.42 
1.22 
4.18 
1.01 
3.99 
1.40 
0.49 
0.61 
6.72 
0.15 
0.07 
2.30 
0.64 
0.42 
9.22 
3.48 
7.13 
6.36 
1.50 
0.88 
1.60 
115.58 
3.74 
2.50 
2.87 
1.15 
0.71 
4.84 
0.96 
0.91 
2.93 
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ID^* 


x^ 


yi. 


R.A. 


Dccl. 


i'775 








ID 


B05 


(pixel) 


(pixel) 


(J2000.0) 


(J2000.0) 


(AB mag) 


Zspec" Ref.<* Zt" 


^mod 




c 


" 


1209 


1038 


13 11 24.979 


-1 20 41.84 


24.451 ± 0.017 


2.6 


4 2.69 ±0.36 


7.18 


18 


a 18.2 


2499 


3112 


13 11 33.066 


-1 20 27.47 


24.774 ±0.024 


1.8 


4 0.16±»:?t 


1.86 




b 18.1 


2130 


1040 


13 11 26.289 


-1 20 00.19 


24.711 ± 0.030 


(1.8) 


•• omtlfl 


4.06 




c 


" 


2029 


2106 


13 11 29.364 


-1 20 27.39 


25.437 ±0.050 


(1.8) 




omtlfi 


12.28 


19 


a 




2404 


2682 


13 11 31.634 


-1 20 22.64 


24.423 ± 0.026 


(2.6) 




0-201^^2 


4.43 




b 


" 


1637 


924 


13 11 25.241 


-1 20 20.05 


24.666 ± 0.026 


(2.6) 




0.20tl:\l 


0.75 




c 




1780 


3082 


13 11 31.958 


-1 20 59.38 


23.472 ±0.011 


(2.6) 




f, 21+0.12 


1.39 




d 




1825 


3088 


13 11 32.040 


-1 20 57.47 


23.842 ±0.017 


2.6 


4 


0.15±0.11 


4.28 


21 


a 




1907 


2714 


13 11 31.027 


-1 20 45.82 


24.470 ±0.034 






1.76t?-f2 


0.50 




b 




1897 


2650 


13 11 30.819 


-1 20 44.91 


23.667 ±0.015 






n r;o + 0.16 


2.06 




c 




1797 


852 


13 11 25.250 


-1 20 11.28 


25.340 ±0.036 






1.78 ±0.27 


7.52 


22 


a 




2407 


2038 


13 11 29.694 


-1 20 08.84 


23.676 ±0.010 


1.7 


4 n n7+°-3* 


0.63 




b 




2127 


2144 


13 11 29.617 


-1 20 23.76 


24.180 ±0.040 


(1.7) 


1 70 + 0.27 


4.48 




c 




1539 


3348 


13 11 32.420 


-1 21 15.94 


23.190 ±0.009 


(1.7) 




n7+''-36 


0.64 


23 


a 




2364 


2005 


13 11 29.533 


-1 20 10.08 


24.640 ±0.023 






Q 09+0.51 

"■'^^-0.09 


1.34 




b 




2134 


2121 


13 11 29.558 


-1 20 22.96 


24.888 ± 0.048 






44+l-''8 
U-44_Q 32 


1.52 




c 




1582 


3408 


13 11 32.662 


-1 21 15.26 


24.537 ±0.018 






0.58tJ:?^ 


1.31 


24 


a 




1485 


2304 


13 11 29.192 


-1 20 56.22 


25.065 ± 0.048 


(2.6) 




O.Sitlil 


0.91 




b 




3039 


2525 


13 11 32.059 


-1 19 50.56 


23.810 ±0.020 


(2.6) 




3.63 ±0.45 


2.53 




c 




3110 


1906 


13 11 30.290 


-1 19 34.21 


24.083 ±0.019 


2.6 


4 


n 99+2.86 


1.22 




d 




2719 


3225 


13 11 33.719 


-1 20 19.91 


24.417 ±0.024 


(2 6) 




2.94 ±0.39 


0.07 


26 


a 




1396 


1008 


13 11 25.153 


-1 20 32.74 


24.880 ±0.019 








1.31 ±0.23 


0.09 




b 




2317 


2621 


13 11 31.326 


-1 20 25.28 


25.334 ±0.032 








0.92 ±0.19 


3.45 




c 




2046 


2389 


13 11 30.242 


-1 20 32.63 


26.525 ±0.200 








qq+2.70 


8.06 


27 


a 




1391 


1017 


13 11 25.173 


-1 20 33.16 


24.976 ±0.019 








-1 n8+0-20 
^■^°-0.90 


1.17 




b 




2334 


2627 


13 11 31.369 


-1 20 24.64 


24.998 ±0.034 








1 11+0-21 


2.58 




c 




2034 


2378 


13 11 30.192 


-1 20 32.94 


27.038 ± 0.225 








1 12+2-°'^ 
-'-■-'-^-0.85 


1.00 


28 


a 




2193 


1677 


13 11 28.301 


-1 20 10.86 


26.717 ±0.090 








c 10+O.6O 
0-l'J_4 36 


1.22 




b 




2062 


3713 


13 11 34.262 


-1 21 00.01 


25.424 ±0.036 








0.72 ±0.17 


6.83 




c 




2034 


2080 


13 11 29.293 


-1 20 26.61 


25.610 ±0.069 








n sn+O-40 

U-SU_o.28 


0.69 


29 


a 




1456 


2329 


13 11 29.226 


-1 20 58.06 


25.805 ± 0.062 


(2.5) 




2.6lt»-f7 


1.75 




b 




3076 


1833 


13 11 30.022 


-1 19 34.20 


24.874 ±0.035 


(2.5) 




3.58 ±0.45 


0.67 




c 




3012 


2577 


13 11 32.177 


-1 19 52.89 


25.169 ±0.073 


(2.5) 




3.61 ±0.45 


1.56 




d 




2694 


3206 


13 11 33.626 


-1 20 20.64 


25.393 ±0.041 


2.5 


4 


3.21 ±0.41 


0.85 


30 


a 




3642 


2362 


13 11 32.420 


-1 19 19.80 


25.938 ±0.055 


(3.0) 




3 Of;+0.46 
''■"'°-3.00 


1.12 




b 




3626 


2623 


13 11 33.185 


-1 19 26.07 


26.032 ± 0.050 


(3.0) 




90+0.42 
■^■^*-3.09 


1.10 




c 


" 


3566 


2809 


13 11 33.662 


-1 19 32.73 


25.814 ±0.045 


3.0 


A 3.0512:2! 


2.37 


31 


a 12.1 


1890 


2333 


13 11 29.852 


-1 20 38.50 


23.985 ± 0.024 


1.83 


2,4 0.70to'il 


0.71 




b HL 


2219 


3312 


13 11 33.274 


-1 20 44.39 


25.373 ±0.037 


(1.83) 


2.oot'l:ll 


1.94 




c 12.4 


2245 


1064 


13 11 26.524 


-1 19 55.49 


25.580 ±0.053 


(1.83) 


0.23+^-^^ 


1.16 




d HL 


1760 


1282 


13 11 26.496 


-1 20 22.08 


25.063 ± 0.030 


(1.83) 


r, oc+o:33 

Z.OO_2 gg 


4.51 


32 


a 


L 


2821 


2671 


13 11 32.191 


-1 20 03.53 


24.083 ±0.013 


(3.0) 


0.15±0.11 


1.97 




b 


L 


2637 


3097 


13 11 33.216 


-1 20 20.90 


25.562 ±0.036 


3.0 


4 3.7lir4? 


1.92 




c 


L 


1416 


2468 


13 11 29.589 


-1 21 02.83 


26.145 ±0.066 


(3.0) 


07+0.43 


0.47 




d 


L 


1796 


2361 


13 11 29.804 


-1 20 43.35 


25.126 ±0.054 


(3.0) 


48+^-^1 


5.90 


33 


a 


L 


1310 


2140 


13 11 28.449 


-1 21 00.66 


27.552 ±0.106 


4.58 


4 0.61+f,-^^ 


2.62 




b 25.2 


2745 


2621 


13 11 31.932 


-1 20 05.91 


26.397 ±0.046 


(4.58) 


3.79 ±0.47 


0.38 


35 


a 


L 


1348 


2159 


13 11 28.560 


-1 20 59.35 


24.632 ± 0.029 


1.9 


4 2.95 ±0.39 


1.91 




b 


L 


2521 


3397 


13 11 33.958 


-1 20 32.52 


24.841 ± 0.031 


(1.9) 


03+301 

^■'^■^-0.03 


9.35 




c 


L 


1906 


2186 


13 11 29.431 


-1 20 34.66 


25.685 ±0.105 


(1.9) 


-1 ^9+0.98 


20.65 


36 


a 


L 


3060 


2352 


13 11 31.566 


-1 19 45.94 


25.639 ± 0.045 


3.0 


A 1 n4+''-98 


0.75 




b 


L 


3049 


2397 


13 11 31.686 


-1 19 47.39 


25.997 ±0.057 


(3.0) 


1 SA+O-SO 
l-84_o.74 


0.46 


40 


a 


L 


2421 


2219 


13 11 30.260 


-1 20 12.04 


27.034 ±0.118 


2.52 


4 l-78tl'58 


0.68 




b 


L 


973 


1545 


13 11 26.176 


-1 21 03.29 


25.876 ±0.050 


(2.52) 


2.88 ±0.38 


1.26 


41 


a 




1453 


1882 


13 11 27.872 


-1 20 48.71 


23.062 ±0.010 








0.55_o:47 


1.57 




b 




2777 


3795 


13 11 35.522 


-1 20 29.38 


25.143 ±0.033 








3.12 ±0.40 


0.66 




c 




1856 


2054 


13 11 28.962 


-1 20 34.12 


99.000 ±28.718 








1 00+3. 83 
^■"'''-0.95 


1.48 


42 


a 




2745 


1541 


13 11 28.672 


-1 19 42.98 


22.969 ±0.015 








0.53_o'23 


0.43 




b 




2885 


2336 


13 11 31.270 


-1 19 53.52 


24.769 ±0.021 








2.68 ±0.36 


1.28 




c 




2403 


3304 


13 11 33.510 


-1 20 35.89 


25.657 ±0.029 








2.81 ±0.39 


2.80 




d 




1261 


2357 


13 11 29.035 


-1 21 07.49 


25.556 ±0.034 








9 c:9+0.35 


0.56 


44 


a 




2030 


1825 


13 11 28.517 


-1 20 21.38 


24.487 ±0.034 








9 r;o + 0.36 
2-5^-0.55 


1.40 
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25 



ID 



45 
46 

48 

49 

50 



ID=' 
B05 



25 



(pixel) 

2060 
2522 
2769 
1964 
1711 
2144 
1660 
2216 
1852 
2145 
1485 
2570 



(pixel) 

3803 
3069 
3852 
2900 

796 
2779 

868 
1849 
3582 
3117 
2910 
2613 



R.A. 
(J2000.0) 



13 11 
13 11 
13 11 
13 11 
13 11 
13 11 
13 11 
13 11 
13 11 
13 11 
13 11 
13 11 



34.531 
32.969 
35.682 
31.669 
24.959 
31.558 
25.104 
28.853 
33.569 
32.580 
31.021 
31.660 



Decl. 
(J2000.0) 

-1 21 02.01 
-1 20 25.51 
-1 20 30.95 
-1 20 47.19 
-1 20 13.98 
-1 20 36.47 
-1 20 17.82 
-1 20 13.47 
-1 21 06.73 
-1 20 43.60 
-1 21 09.08 
-1 20 13.66 



« 775 
(AB mag) 



Ref. 



25.308 ± 
25.138 ± 
26.066 ± 
25.121 ± 
25.674 ± 
25.215 ± 
27.085 ± 
25.602 ± 
26.285 ± 
26.691 ± 
26.595 ± 
28.007 ± 



0.042 
0.058 
0.041 
0.040 
0.053 
0.047 
0.111 
0.049 
0.060 
0.061 
0.047 
0.182 



7+0. 

3.70 ± 
0.62 ± 


2 
0, 
1 
1 

+0, 
2 

+2, 
0, 

+0, 
2 

+0. 
0, 

+0, 
0, 

1 

2 



2.27^ 



0.76 
0.20 
1.41 
2.47 
0.2O 
2.60 
4.08 
3.25 
2.23 



.11 

0.46 

0.16 

.37 
59 
50 
12 
19 
27 
66 
10 
73 
12 
40 
32 
50 
52 
99 
42 
74 
02 



■^moa 

0.51 
1.23 
5.33 
135.43 
1.50 
0.23 
1.25 
0.65 
0.90 
0.42 
0.29 
1.71 



and 
b 



Quotation marks indicate agreement with the B05 numbering system (where our la is their 1.1, etc.). ''L" indicates a system added by LOG, 

" a system added in this paper. "HL" indicates systems that were adde d by H06, but for which we use L06's numbering. 
Pixel coordinates (0. 05/pix) in the original APSlS-reduced ACS images (see E3.5i. 
Values in parentheses a re assumed from other multip le imagefs) . 



1) 



Frye et al. 



I 2002 I ; 2) B05; 3) 



Frye et al. 



I 2007 I ; 4) LOT 



Reference for Zs 

Bayesian photometric redshift and 95% confidence interval. Based on fits to photometry across 12 filters without subtraction of cluster galaxies. 



Poorness of BPZ fit: observed vs. model fl 



uxes wr 



th uncertainties 



TABLE 6 
Extra Knots Identified in the Multiple Images 



ID^ 



(pixel) 



y" 

(pixel) 



R.A. 
(J2000.0) 



ID numbers as in Table Isl with the third column giving the "knot number'' 

Pixel coordinates (0''05/pix) in the original APSIS-rcduced ACS images (sec 3.5 i. 



Decl. 
(J2000.0) 



4 


a 


1 


1818 


3116 


13 11 32.114 


-1 20 58.38 




b 


1 


1395 


2832 


13 11 30.658 


-1 21 11.50 




c 


1 


2545 


2328 


13 11 30.765 


-1 20 08.74 




d 


1 


1485 


1336 


13 11 26.269 


-1 20 35.67 


6 


a 


1 


3093 


2082 


13 11 30.798 


-1 19 38.71 




a 


2 


3092 


2075 


13 11 30.775 


-1 19 38.61 




b 


1 


2803 


3067 


13 11 33.361 


-1 20 12.75 




b 


2 


2801 


3058 


13 11 33.331 


-1 20 12.65 


8 


a 


1 


1956 


3098 


13 11 32.255 


-1 20 51.75 




a 


2 


1996 


3108 


13 11 32.342 


-1 20 50.15 




b 


1 


1628 


3000 


13 11 31.495 


-1 21 04.52 




b 


2 


1560 


2961 


13 11 31.281 


-1 21 06.77 


10 


a 


1 


2179 


3564 


13 11 33.978 


-1 20 51.55 




b 


1 


2134 


1634 


13 11 28.088 


-1 20 12.62 


16 


a 


1 


1897 


1732 


13 11 28.048 


-1 20 25.43 




a 


2 


1864 


1712 


13 11 27.941 


-1 20 26.50 




b 


1 


1943 


1958 


13 11 28.795 


-1 20 28.14 




b 


2 


1946 


2005 


13 11 28.942 


-1 20 29.00 




c 


1 


2327 


3643 


13 11 34.426 


-1 20 46.53 


18 


a 


1 


2095 


3555 


13 11 33.832 


-1 20 55.16 




b 


1 


2206 


1663 


13 11 28.277 


-1 20 09.97 


19 


a 


1 


2384 


2716 


13 11 31.708 


-1 20 24.27 




a 


2 


2421 


2644 


13 11 31.543 


-1 20 21.06 




b 


1 


1657 


910 


13 11 25.227 


-1 20 18.84 




b 


2 


1613 


941 


13 11 25.258 


-1 20 21.49 




d 


1 


1834 


3090 


13 11 32.059 


-1 20 57.11 


31 


a 


2 


2807 


2022 


13 11 30.212 


-1 19 50.39 




a 


3 


2814 


2070 


13 11 30.366 


-1 19 51.09 




b 


3 


2233 


3311 


13 11 33.290 


-1 20 43.73 




c 


2 


1199 


2349 


13 11 28.923 


-1 21 10.12 




c 


3 


1207 


2358 


13 11 28.961 


-1 21 09.95 




d 


2 


1762 


1288 


13 11 26.517 


-1 20 22.11 




d 


3 


1785 


1276 


13 11 26.513 


-1 20 20.82 



26 
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TABLE 7 
Multiple Image Systems 





Number 


Knots 


Input Spectroscopic 


Input 


Output 


ID 


of images 


per image 


order redshift Ref.'' redshift'^ 


redshift '^ 


1 


6 




1 3.04 


2,3 3.04 


3.04 


2 


5 




2 2.53 


3 2.54 


2.54 


3 


3 




23 




5.47 ±0.7 


5.490 


4 


5 




7 


L.l 


4 1.1 


1.1 


5 


3 




13 


2.6 


4 2.6 


2.6 


6 


4 




12 


L.l 


4 1.1 


1.1 


7 


2 




5 4.87 


1 4.87 


4.87 


8 


4 




10 




1.8 ±0.5 


2.224 


9 


4 




4 




5.16 ±0.5 


4.707 


10 


3 




18 1.83 


4 1.8 


1.8 


11 


3 




25 


2.5 


4 2.5 


2.5 


12 


6 




30 1.83 


2,4 1.83 


1.83 


13 


3 




31 




1.02 ±0.5 


1.532 


14 


2 




32 


3.4 


4 3.4 


3.4 


15 


3 




19 


L.8 


4 (10) 


1.8 


16 


3 




14 




2.01 ±0.5 


1.635 


17 


3 




26 


2.6 


4 2.6 


2.6 


18 


3 




20 


L.8 


4 (10) 


1.8 


19 


4 




11 


2.6 


4 2.6 


2.6 


21 


3 




15 




1.78 ±0.5 


1.065 


22 


3 




16 


L.7 


4 1.7 


1.7 


23 


3 




17 




■■ (22) 


1.7 


24 


4 




8 


2.6 


4 2.6 


2.6 


26 


3 




21 




■■ (27) 


2.388 


27 


3 




22 




1.74 ±0.5 


2.388 


28 


3 




24 




5.45 ±0.7 


4.451 


29 


4 




9 


2.5 


4 2.5 


2.5 


30 


3 




33 


i.O 


4 3.0 


3.0 


31 


4 




3 1.83 


2,4 1.8 


1.8 


32 


4 




6 


i.O 


4 3.0 


3.0 


33 


2 




27 4.58 


4 4.5 


4.5 


35 


3 




28 


L.9 


4 1.9 


1.9 


36 


2 




34 


3.0 


4 3.0 


3.0 


40 


2 




29 2.52 


4 2.5 


2.5 


41 


3 




37 




2.5 ±2.0 


5.357 


42 


4 




35 






2.0 ±2.0 


2.222 


44 


2 




38 






2.0 ±2.0 


3.210 


45 


2 




39 






2.5 ±2.0 


5.359 


46 


2 




40 






2.5 ±2.0 


1.714 


48 


2 




41 






2.0 ±2.0 


1.456 


49 


2 




42 






2.0 ±2.0 


1.645 


50 


3 




36 




2.5 ±2.0 


3.029 



Note. — Our mass model was built by iteratively adding our multiple image systems in the order given here. Any additional knots were added 
immediately after inclusion of the "main knot" for that system. 

'^ Reference for Zsp^^: l) |Frye et al]|2002| ; 2) BOS; 3) |Frye et al] l |2007| ; 4) L07 

Input rcdshifts with rough uncertainties where allowed. (Spectroscopic redshifts arc given zero uncertainty.) Rcdshifts of some systems are tied 
to other systems with numbers given in parentheses. (For example, the redshift of system 26 is tied to that of system 27 in the optimization.) 

Final model redshift, optimized if input value included uncertainty. 



APPENDIX 
A. DARK MATTER HALO MASS PROFILES 



Clus ters formed in simula tions have m a.ss profiles generally well described by NFW ( Navarro et al.|19"96 l or Einasto / 
Scrsic ("Navarro et al. 12004') profi les (e.g., Navarro ct al.'201 ()|). The latter were found to yield superior hts to simulated 



halos (Navarro ct al 
in 



2004; Mcrritt et al. 



Coel (|2010|) 
NT 



2005, 2006; Navarro et al.|[20l0 ) . Details on these fitting forms can be found 



'he INi^'W profile gives mass density as a function of radius in three dimensions: 



p{r) = 



(r/r,)(l + r/r,)2 



(Al) 



The logarithmic slope increases from -1 at the core to -3 in the outskirts. The NFW concentration parameter is given 
as c = ryir/Vs- By definition, the average density within the virial radius r^^^ is Aj, ss 115 ti mes Prrit = SH'^ ( z) / (SttG ) 
for a collapsed virialized sphere at z = 0.186, as evaluated according to formulae given by Nakamura fc Suto| ( |1997 ). 
For historical reasons, Ac « 200 is often used, so we report concentration values for both overdensities. C onversion 
between the two is given by C115 w 1.234c2oo + 0.172. Lensing properties of the NFW profile can be found in Golse & 
|Kneibj(;2002j). 
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Fig. 16. — Multiple images in Abell 1689. For each image, identified knots are labeled with numbers and crosses. Green lines are 
isophotal apertures, altered where necessary. Yellow boxes outlining the IDs indicate color image stamps that were produced from our 
galaxy-subtracted images. All other color stamps are taken from the STScI g'r' z' color image. Coordinates in pixels are given along the 
axes. Each pixel measures 0''05 across. 



10a ^^ - 10b J^ - 10^^^- 



16a 



^ I "A 16b %^ I ="'| 16c il 



Fig. 17. — Multiple images delensed to the source plane. Knots are labeled as in Fig. |16| but now constrained knots are colored green 
and unconstrained knots are colored pink. Within each system, all of the constrained knots align exactly among the delensed images. 
Unconstrained knots are generally well aligned but not always. Panels are plotted to the same scale within each system. Yellow boxes 
outlining the IDs indicate color image stamps that were produced from our galaxy-subtracted images. All other color stamps are taken 
from the STScI g'r' z' color image. Coordinates in pixels are given along the axes. Each pixel measures O'.'OS across. 



The Sersic ( 1968 1 profile has a co ntinuously varying s lope t hat asymptotes toward zero a t very small radius. Its 



use as a mass protiie is descri bed by 



calculated by Cardone (2004) and E' 



Sersic profile of projected mass is given as: 



E(i?) 



Graham fc Driver (2005) and Terzic & Graham (2005) with lensing pro perties 
iasdottir & Moilerl ( 2007 1 . These details are summarized in Coe (2010). Our 



, exp' 



l/r. 



(A2) 



There are three free pa rameters: Sg, i?e, and n, with 6„ being a function of n. This profile was found (Navarro 
et al.||20d4 Merritt et al. 2005, 2006 ) to yield better fits to a wide range of simulated dark matter halos than did the 
generalized IN l^'W protiie ( |Zhao|1996[|Wyithe et al.|2001 ), which has an equal number (3) of free parameters, including 
the central slope. 

A distinction has been made between the Sersic profile described above and the Einasto profile. The latter has a 
similar form but is defined as a function of 3-D density p(r) rather than projected surface density S(i?) — K( R)T,crit- 
Performanc e of the two forms is similar with the Sersic profile providing a slightly be tter fit to cluster halos ( [Merritt 
et al.|[2006). These and other mass profiles are described in more detail in Goe (2010). 



B. MULTIPLE IMAGES 

In this section we display color images of the multiple images used in this work. Fig. [16] displays the images as they 
appear and Fig. 17 displays the same images delensed, as they would appear if Abell 1689 did not exist (according 
to our best mass model). The 55 constrained positions are labeled in green in Fig. 17 These delensed positions 



match perfectly among each set of multiple images. We have identified other knots, labeled in pink, which we have 
not constrained. We have doubts about these identifications (including a case or two in which it appears that we were 
obviously mistaken). 

Examples of the image stamps are given here. The remaining images (for all 42 image systems) can be found in the 
online version. 



